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We analyze and compare candidate crystal structures for the crystalline color superconducting 
phase that may arise in cold, dense but not asymptotically dense, three-flavor quark matter. We 
determine the gap parameter A and free energy Q{A) for many possible crystal structures within a 
Ginzburg-Landau approximation, evaluating Q{A) to order A^. In contrast to the two-flavor case, 
we find a positive A® term and hence an S7(A) that is bounded from below for all the structures 
that we analyze. This means that we are able to evaluate A and £7 as a function of the splitting 
between Fermi surfaces for all the structures we consider. We find two structures with particularly 
robust values of A and the condensation energy, within a factor of two of those for the CFL phase 
which is known to characterize QCD at asymptotically large densities. The robustness of these 
phases results in their being favored over wide ranges of density. However, it also implies that 
the Ginzburg-Landau approximation is not quantitatively reliable. We develop qualitative insights 
into what makes a crystal structure favorable, and use these to winnow the possibilities. The two 
structures that we find to be most favorable are both built from condensates with face-centered 
cubic symmetry: in one case, the (ltd) and (us) condensates are separately face centered cubic; in 
the other case (ud) and (us) combined make up a face centered cube. 

PACS numbers: 12.38.-t, 26.60.+C, 12.38.Mh, 74.20.-z 



I. INTRODUCTION 



Quantum chromodynamics predicts that at densi- 
ties that are high enough that baryons are crushed 
into quark matter, the quark matter that results fea- 
tures pairing between quarks at low enough temper- 
atures, meaning that it is in one of a family of pos- 
sible color superconducting phases 0. The essence 
of color superconductivity is quark pairing driven 
by the BCS mechanism, which operates whenever 
there are attractive interactions between fermions at 
a Fermi surface • The interaction between quarks 
in QCD is strong and is attractive between quarks 
that are antisymmetric in color, so we expect cold 
dense quark matter to exhibit color superconductiv- 
ity. If color superconducting quark matter occurs in 
nature, it lies within compact stars. Except during 
the first few seconds after their birth in supernovae, 
these stars have temperatures well below the tens of 
MeV. This implies that if these stars feature quark 
matter cores, these cores will be color superconduc- 
tors, and justifies us in restricting our investigation 
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to T = throughout this paper. 

Wc shall only consider Cooper pairs whose pair 
wave function is antisymmetric in Dirac indices 
— the relativistic generalization of zero total spin. 
(Other possibilities have been investigated 0, 13, S 
0, Q and found to be less favorable.) This in turn 
requires antisymmetry in flavor, meaning in partic- 
ular that the two quarks in a Cooper pair must have 
different flavor. 

It is by now well-established that at sufficiently 
high densities, where the up, down and strange 
quarks can be treated on an equal footing and the 
disruptive effects of the strange quark mass can be 
neglected, quark matter is in the color-flavor locked 
(CFL) phase, in which quarks of all three colors 
and all three flavors form conventional Cooper pairs 
with zero total momentum, and all fermionic exci- 
tations are gapped, with the gap parameter Ag ~ 
10—100 MeV 0,0 ■ However, even at the very center 
of a compact star the quark number chemical poten- 
tial n cannot be much larger than 500 MeV, meaning 
that the strange quark mass Afs (which is density de- 
pendent, lying somewhere between its vacuum cur- 
rent mass of about 100 MeV and constituent mass of 
about 500 MeV) cannot be neglected. Furthermore, 
bulk matter, as relevant for a compact star, must 
be in weak equilibrium and must be electrically and 
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color neutral H, S H IS ^3 ■ AH these factors work 
to separate the Fermi momenta of the three different 
flavors of quarks, and thus disfavor the cross-species 
BCS pairing that characterizes the CFL phase. If 
we imagine beginning at asymptotically high densi- 
ties and reducing the density, and suppose that CFL 
pairing is disrupted by the heaviness of the strange 
quark before color superconducting quark matter is 
superseded by baryonic matter, the CFL phase must 
be replaced by some phase of quark matter in which 
there is less, and less symmetric, pairing. 

Within a spatially homogeneous ansatz, the next 
phase down in density is the gapless CFL (gCFL) 
phase [IllllllSIIlllllIlllllIi- In this phase, 
quarks of all three colors and all three flavors still 
form ordinary Cooper pairs, with each pair hav- 
ing zero total momentum, but there are regions of 
momentum space in which certain quarks do not 
succeed in pairing, and these regions are bounded 
by momenta at which certain fermionic quasipar- 
ticles are gapless. This variation on BCS pairing 
— in which the same species of fermions that pair 
feature gapless quasiparticles — was first proposed 
for two flavor quark matter [T9| and in an atomic 
physics context [23|. In all these contexts, how- 
ever, the gapless paired state turns out in general 
to suffer from a "magnetic instability" : it can lower 
its energy by the formation of counter-propagating 
currents |^ [2^ In the atomic physics con- 

text, the resolution of the instability is phase sep- 
aration, into macroscopic regions of two phases in 
one of which standard BCS pairing occurs and in 
the other of which no pairing occurs 0, [2^ |2^ . In 
three-flavor quark matter, where the instability of 
the gCFL phase has been established in Refs. p2 |. 
phase coexistence would require coexisting compo- 
nents with opposite color charges, in addition to 
opposite electric charges, making it very unlikely 
that a phase separated solution can have lower en- 
ergy than the gCFL phase 0, 0|. Furthermore, 
color superconducting phases which are less sym- 
metric than the CFL phase but still involve only 
conventional BCS pairing, for example the much- 
studied 2SC phase in which only two colors of 
up and down quarks pair 0, l27l ^3 but includ- 
ing also many other possibilities |29j, cannot be 
the resolution of the gCFL instability It 
seems likely, therefore, that a ground state with 
counter-propagating currents is required. This could 
take the form of a crystalline color superconduc- 
tor [iSlllllllMllllllllllllllllsSliilil 



— the QCD analogue of a form of non-BCS pair- 
ing first considered by Larkin, Ovchinnikov, Fulde 
and Ferrell given that the CFL phase itself 

is likely augmented by kaon condensation 0i Q| > 
it could take the form of a phase in which a CFL 
kaon condensate carries a current in one direction 
balanced by a counter-propagating current in the 
opposite direction carried by gapless quark quasi- 
particles |45l |4^ . This meson supercurrent phase 
has been shown to have a lower free energy than the 
gCFL phase. 

Our purpose in this paper is to analyze and com- 
pare candidate crystal structures for three-flavor 
crystalline color superconductivity. The investiga- 
tion of crystalline color superconductivity in three- 
flavor QCD was initiated in Ref. [s^. Although 
such phases seem to be free from magnetic insta- 
bility |40| . it remains to be seen whether such a 
phase can have a lower free energy than the me- 
son current phase, making it a possible resolution to 
the gCFL instability. The simplest "crystal" struc- 
tures do not suffice [sol |4]] | , but experience in the 
two-flavor context |3a| suggests that realistic crystal 
structures constructed from more plane waves will 
prove to be qualitatively more robust. Our results 
confirm this expectation. 

Determining the favored crystal structure (s) in 
the crystalline color superconducting phase(s) of 
three-flavor QCD requires determining the gaps and 
comparing the free energies for very many candidate 
structures, as there are even more possibilities than 
the many that were investigated in the two-flavor 
context [33. As there, we shall make a Ginzburg- 
Landau approximation. This approximation is con- 
trolled if A ^ Aq, where A is the gap parameter 
of the crystalline color superconducting phase itself 
and Ao is the gap parameter in the CFL phase that 
would occur if Ms were zero. We shall find that 
the most favored crystal structures can have A/Aq 
as large as ~ 1/2, meaning that we are pushing the 
approximation hard and so should not trust itquan- 
titatively. In earlier work with Mannarelli |41|, we 
analyzed a particularly simple one parameter family 
of "crystal" structures in three-flavor quark matter, 
simple enough that we were able to do do the anal- 
ysis both with and withovit the Ginzburg-Landau 
approximation. We found that the approximation 
works when it should and that, at least for the simple 
crystal structures we analyzed in Ref. when it 
breaks down it always underestimates the gap A and 
the condensation energy. Furthermore, we found 
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that the Ginzburg-Landau approximation correctly 
determines which crystal structure among the one 
parameter family that we analyzed in Ref. |4ll | has 
the largest gap and lowest free energy. 

We shall work throughout in a Nambu-Jona- 
Lasinio (NJL) model in which the QCD interaction 
between quarks is replaced by a point-like four-quark 
interaction, with the quantum numbers of singlc- 
gluon exchange, analyzed in mean field theory. This 
is not a controlled approximation. However, it suf- 
fices for our purposes: because this model has at- 
traction in the same channels as in QCD, its high 
density phase is the CFL phase; and, the Fermi sur- 
face splitting effects whose qualitative consequences 
we wish to study can all be built into the model. 
Note that we shall assume throughout that Aq <C fi. 
This weak coupling assumption means that the pair- 
ing is dominated by modes near the Fermi surfaces. 
Quantitatively, this means that results for the gaps 
and condensation energies of candidate crystalline 
phases are independent of the cutoff in the NJL 
model when expressed in terms of the CFL gap Aq: 
if the cutoff is changed with the NJL coupling con- 
stant adjusted so that Aq stays fixed, the gaps and 
condensation energies for the candidate crystalline 
phases also stay fixed. This makes the NJL model 
valuable for making the comparisons that are our 
goal. 

We shall consider crystal structures in which there 
arc two condensates 

{ud) - A3 ^ exp (2^q;^ • r) 

a 

(us) ^ A2 ^ exp (2^q^ r) . (1) 

a 

As in Refs. [S^^J, and as we explain in Section II, 
we neglect (ds) pairing because the d and s Fermi 
surfaces are twice as far apart from each other as 
each is from the intervening u Fermi surface. Were 
we to set A2 to zero, treating only (ud) pairing, we 
would recover the two-flavor Ginzburg-Landau anal- 
ysis of Ref. ■ There, it was found that the best 
choice of crystal structure was one in which pairing 
occurs for a set of eight q|'s pointing at the corners 
of a cube in momentum space, yielding a condensate 
with face-centered cubic symmetry. The analyses 
of three-flavor crystalline color superconductivity in 
Refs. 11^01 introduce nonzero A2, but made the 
simplifying ansatz that pairing occurs only for a sin- 
gle qa and a single q2 . Wc consider crystal structures 
with up to eight qg's and up to eight q2's. 



We shall evaluate the free energy ri(A2,A3) for 
each crystal structure, in a Ginzburg-Landau expan- 
sion in powers of the A's. We work up to order Aj A' 
with p+q = 6. At sextic order, we find that fl{A, A) 
is positive for large A for all the crystal structures 
that we investigate. This is in marked contrast to 
the results of Ref. [sH , which showed that many two- 
flavor crystal structures have negative sextic terms, 
with free energies that are unbounded from below 
when the Ginzburg-Landau expansion is stopped at 
sextic order. Because we find positive sextic terms, 
we are able to use our sextic Ginzburg-Landau ex- 
pansion to evaluate A and il(A, A) for all the struc- 
tures that we analyze. 

The two crystal structures that we argue are most 
favorable are both related to the face-centered cube 
of Ref. [33, but in different ways. In the first, which 
we denote "CubeX" in Section VI, there are four 
qf 's and four qj's which together point at the eight 
corners of a cube in momentum space. In the second, 
denoted "2Cube45z" in Section VI, there are eight 
qg's and eight qj's which each point at the eight cor- 
ners of a cube in momentum space, the two cubes 
rotated relative to each other by 45 degrees about 
an axis perpendicular to their faces. To a large de- 
gree, our argument that these two structures are the 
most favorable relies only on two qualitative inputs. 
First, if either the set of {q2}'s or the set of {qgj's 
yields a (us) or a (ud) condensate whose free energy, 
viewed in isolation as a two-flavor problem and eval- 
uated as in Ref. [s^, is unfavorable, then the three- 
flavor condensate is unfavorable. Thus, we can use 
all the qualitative results of Ref. [s^. Second, the 
free energy of a candidate three-flavor crystal struc- 
ture becomes less favorable the closer any comes 
to the antipodes of any qg. This second result is 
foreshadowed in the results of Refs. H^^J, and the 
results of Ref. indicate that it is valid beyond 
the Ginzburg-Landau approximation. We shall see 
in Section VI that these two qualitative lessons are 
sufficient to winnow the space of candidate crystal 
structures down to the two that our calculational re- 
sults, also described in Section VI, demonstrate are 
indeed the most favorable. 

We find that several of the crystal structures that 
we consider have gap parameters that can be as 
large as Ao/3, and that one of them (the CubeX 
structure) has A/Aq that reaches 1/2. The robust- 
ness of these crystalline condensates thus pushes the 
Ginzburg-Landau approximation that we have used 
in the derivation of our results to the edge of its 
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regime of quantitative reliability. As we discussed 
above, the analysis of Ref. ^J] shows that for simpler 
crystal structures qualitative results obtained within 
this approximation remain valid when the approxi- 
mation has broken down quantitatively. We expect 
this to be so also for the more realistic, and compli- 
cated, crystal structures that we have constructed, 
but a demonstration would require their analysis 
without making a Ginzburg-Landau approximation, 
something we do not attempt here. 

We find that the two crystal structures which we 
argue are most favorable have large condensation en- 
ergies, easily 1/3 to 1/2 of that in the CFL phase 
with Ms = 0, which is SAq^^/tt^. This is remark- 
able, given the only quarks that pair are those lying 
on (admittedly many) rings on the Fermi surfaces, 
whereas in the CFL phase with = pairing oc- 
curs over the entire u, d and s Fermi surfaces. 

The gapless CFL (gCFL) phase provides a use- 
ful comparison at nonzero Mg. For 2Ao < M"^ / ji < 
5.2Ao, model analyses that are restricted to isotropic 
phases predict a gCFL phase 0, 0, finding 
this phase to have lower free energy than either the 
CFL phase or unpaired quark matter. However, 
this phase is unstable to the formation of current- 
carrying condensates [U, [S^, IH, ES IH and so it 
cannot be the ground state. The true ground state 
must have lower free energy than that of the gCFL 
phase, and for this reason the gCFL free energy pro- 
vides a useful benchmark. We find that three-flavor 
crystalline color superconducting quark matter has 
a lower free energy than both gCFL quark matter 
and unpaired quark matter within a wide regime of 
density. For 

2.9Ao < ^ < 10.4Ao (2) 
M 

the crystalline phase with one or other of the two 
crystal structures that we argue are most favorable 
has lower free energy (greater condensation energy) 
than CFL quark matter, gCFL quark matter, and 
unpaired quark matter. (See Fig. |7|in Section VI.) 
This window in parameter space is in no sense nar- 
row. Our results therefore indicate that three-flavor 
crystalline quark matter will occur over a wide range 
of densities, unless, that is, the pairing between 
quarks is so strong (that is, Aq is so large making 
M^/Ag so small) that quark matter is in the CFL 
phase all the way down to the density at which quark 
matter is superseded by nuclear matter. 

However, our results also indicate that unless the 



Ginzburg-Landau approximation is underestimating 
the condensation energy of the crystalline phase by 
about a factor of two, there is a fraction of the 
"gCFL window" (with 2Ao < M^/fi < 2.9Ao, in the 
Ginzburg-Landau approximation) in which no crys- 
talline phase has lower free energy than the gCFL 
phase. This is thus the most likely regime in which 
to find the current-carrying meson condensates of 
Refs. miil. 

Our paper is organized as follows. In Section ITTl 
we shall specify the model we use and the simplifying 
assumptions we make, valid for A ^ Aq. Along the 
way we review relevant aspects of two-flavor color 
superconductivity. We shall also define our ansatz 
for the crystalline condensates more precisely than 
in Eq. . Much of Section ^1 closely follows our 
earlier paper in collaboration with Mannarelli [Zlj . 
One simplifying assumption that we make is that 
A2 and A3 are equal in magnitude, an assumption 
which is related to how electric neutrality is main- 
tained. In Appendix A, we use our results to confirm 
the validity of this assumption. In Section III we in- 
troduce the Ginzburg-Landau expansion of the free 
energy, deferring the derivation of the expressions for 
the Ginzburg-Landau coefficients to Section IV and 
their evaluation to Section V. We give our results in 
Section VI, and discuss their implications for future 
work in Section VII. 



II. MODEL, SIMPLIFICATIONS AND 
ANSATZ 

A. Neutral unpaired three-flavor quark matter 

We shall analyze quark matter containing mass- 
less u and d quarks and s quarks with an effective 
mass Ms- (Although the strange quark mass can be 
determined self-consistently by solving for an (ss) 
condensate 0, 0, 0| , we shall leave this to future 
work and treat Mg as a parameter.) The Lagrangian 
density describing this system in the absence of in- 
teractions is given by 

Co = ^» {^Kf~ + t^tf 70) ^« , (3) 

where i,j — 1,2,3 are flavor indices and a,/3 = 
1,2,3 are color indices and we have suppressed the 
Dirac indices, where M^^^ = (5"'' diag(0, 0, A/s)y is 

the mass matrix, where d^f = dS"^Sij and where 
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the quark chemical potential matrix is given by 



V3 



«/3 



Q/3 

(4) 

with Q = diag(2/3, —1/3, —1/3)^ the quark electric- 
charge matrix and T3 and Tg the GcU-Mann matrices 
in color space. We shall quote results at quark num- 
ber chemical potential ^ = 500 MeV throughout. 

In QCD. /ie, /i3 and /xg are the zeroth components 
of electromagnetic and color gauge fields, and the 
gauge field dynamics ensure that they take on values 
such that the matter is neutral [a,|43i satisfying 



= 



(5) 



with n the free energy density of the system. In the 
NJL model that we shall employ, in which quarks 
interact via four-fermion interactions and there are 
no gauge fields, we introduce /ie, /13 and /ig by hand, 
and choose them to satisfy the neutrality constraints 
(jSJ. The assumption of weak equilibrium is built 
into the calculation via the fact that the only flavor- 
dependent chemical potential is jif. , ensuring for ex- 
ample that the chemical potentials of d and s quarks 
with the same color must be equal. Because the 
strange quarks have greater mass, the equality of 
their chemical potentials implies that the s quarks 
have smaller Fermi momenta than the d quarks in 
the absence of BCS pairing. In the absence of pair- 
ing, then, because weak equilibrium drives the mas- 
sive strange quarks to be less numerous than the 
down quarks, electrical neutrality requires a /ig > 0, 
which makes the up quarks less numerous than the 
down quarks and introduces some electrons into the 
system. In the absence of pairing, color neutrality is 
obtained with /j,3 = /xg = 0. 

The Fermi momenta of the quarks and electrons 
in quark matter that is electrically and color neutral 
and in weak equilibrium are given in the absence of 
pairing by 



Pf 



Pf 

Pf 
Pf 



Me 

3 
3 



3 



-AH 
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Ml 
2^l 



Me 



(6) 



where we have simplified p^p upon assuming that Mg 
and Me are small compared to /i by working only to 



linear order in //e and M^. The free energy of the 
noninteracting quarks and electrons is given by 



?^{plf + i{v%Y + {p%)' 



127r2 



Pf 



p^dp ( \J p^ + Mg — n — 



unpaired 



To this order, electric neutrality requires 

Me = -j— , 

4m 



(7) 



(8) 



yielding 



Pf 



Pf 
Pf 
Pf 



= M 



Ml 
12/i 
Ml 
6m 

12/1 



Pf 



Pf 



Ml 
4m 

Ml 
4m 



Ml 
4/i 



(9) 



We see from (|SJ that to leading order in and 
Me, the effect of the strange quark mass on unpaired 
quark matter is as if instead one reduced the strange 
quark chemical potential by M^/(2/i). We shall 
make this approximation throughout. The correc- 
tions to this approximation in an NJL analysis of 
a two-fiavor crystalline color superconductor have 
been evaluated and found to be small [s^l, and we 
expect the same to be true here. Upon making this 
assumption, we need no longer be careful about the 
distinction between pf's and m'Sj as we can simply 
think of the three flavors of quarks as if they have 
chemical potentials 



Md = Mu + 2(5m3 

M« = Pf 

Ms = M« - 2(5/i2 



with 



(5M3 = 5M2 



Ml 
8m 



= <5m , 



(10) 



(11) 
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where the choice of subscripts indicates that 25^2 is 
the spUtting between the Fermi surfaces for quarks 1 
and 3 and 26^3 is that between the Fermi surfaces for 
quarks 1 and 2, identifying u,d,s with 1,2,3. (The 
prefactor 2 in the equations defining the (5/z's is cho- 
sen to agree with the notation used in the analysis 
of crystaUine color superconductivity in a two flavor 
model Is^l , in which the two Fermi surfaces were de- 
noted hy fi ± Sfi meaning that they were separated 
by 2(5^.) 

Note that the equality of 6^,2 and Sfi^ is only valid 
to leading order in M^; at the next order, /ie = 
M2/(4^) - Mf/{4:8^l^) and 6^3 = Me/2 while 6^12 = 
(5/i3 4- / {16fi^). In Section V, we will utilize the 
fact that Sfi2 and Sfi3 are close to equal, but not 
precisely equal. 



B. BCS pairing and neutrality 

As described in Refs. [iSlillil, BCS pairing 
introduces qualitative changes into the analysis of 
neutrality. For example, in the CFL phase /ig = 
and is nonzero and of order M^//x. This arises 
because the construction of a phase in which BCS 
pairing occurs between fermions whose Fermi sur- 
face would be split in the absence of pairing can 
be described as follows. First, adjust the Fermi sur- 
faces of those fermions that pair to make them equal. 
This costs a free energy price of order Sfj.'^fj,^. And, it 
changes the relation between the chemical potentials 
and the particle numbers, meaning that the /i's re- 
quired for neutrality can change qualitatively as hap- 
pens in the CFL example. Second, pair. This yields 
a free energy benefit of order Aq/x^, where Aq is the 
gap parameter describing the BCS pairing. Hence, 
BCS pairing will only occur if the attraction between 
the fermions is large enough that Aq > Sfj,. In the 
CFL context, in which (ud), (us) and (ds) pairing is 
fighting against the splitting between the d, u and 
s Fermi surfaces described above, it turns out that 
CFL pairing can occur if Ao > 4Sfi = M^/{2n) jnl, 
a criterion that is reduced somewhat by kaon con- 
densation which acts to stabilize CFL pairing |4^ . 

In this paper we are considering quark matter at 
densities that are low enough {fi < M^/(2Ao)) that 
CFL pairing is not possible. The gap parameter Aq 
that would characterize the CFL phase if and 
were zero is nevertheless an important scale in our 
problem, as it quantifies the strength of the attrac- 
tion between quarks. Estimates of the magnitude of 



Aq are typically in the tens of MeV, perhaps as large 
as 100 MeV |l|. We shall treat Ag as a parameter, 
and quote results for Aq = 25 MeV, although as 
we shall show in Section VI. E our results can easily 
be scaled to any value of Aq as long as the weak- 
coupling approximation Aq <C is respected. 

C. Crystalline color superconductivity in 
two-flavor quark matter 

Crystalline color superconductivity can be 
thought of as the answer to the question: "Is there 
a way to pair quarks at differing Fermi surfaces 
without first equalizing their Fermi momenta, given 
that doing so exacts a cost?" The answer is "Yes, 
but it requires Cooper pairs with nonzero total 
momentum." Ordinary BCS pairing pairs quarks 
with momenta p and — p, meaning that if the Fermi 
surfaces are split at most one member of a pair 
can be at its Fermi surface. In the crystalline color 
superconducting phase, pairs with total momentum 
2q condense, meaning that one member of the 
pair has momentum p -I- q and the other has 
momentum — p + q for some p USUI. Suppose 
for a moment that only u and d quarks pair, 
making the analyses of a two-flavor model found in 
Refs. [23, EH IH ElB El El El m (and really 
going back to Ref. [ij) valid. We sketch the results 
of this analysis in this subsection. 

The simplest "crystalline" phase is one in which 
only pairs with a single q condense, yielding a con- 
densate 

{ipu{x)Cj5'ilJd{x)) oc A exp(2iq • r) (12) 

that is modulated in space like a plane wave. (Here 
and throughout, we shall denote by r the spa- 
tial three-vector corresponding to the Lorentz four- 
vector X.) Assuming that A ^ Sfi <^ fi, the en- 
ergetically favored value of |q| = g turns out to be 
q = 1^6 fi, where the proportionality constant rj is 
given by 77 = 1.1997 [S^, ■ If V were 1, then the 
only choice of p for which a Cooper pair with mo- 
menta (— p-l-q, p-l-q) would describe two quarks each 
on their respective Fermi surfaces would correspond 
to a quark on the north pole of one Fermi surface and 
a quark on the south pole of the other. Instead, with 
7] > 1, the quarks on each Fermi surface that can 
pair lie on one ring on each Fermi surface, the rings 
having opening angle ■00 = 2cos~^(l/77) = 67.1°. 
The energetic calculation that determines rj can be 
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thought of as balancing the gain in pairing energy 
as 77 is increased beyond 1, allowing quarks on larger 
rings to pair, against the kinetic energy cost of 
Cooper pairs with greater total momentum. If the 
A/Sfj, — *■ Ginzburg-Landau limit is not assumed, 
the pairing rings change from circular lines on the 
Fermi surfaces into ribbons of thickness ^ A and an- 
gular extent ~ A/Sfi. The condensate (|12ll carries a 
current, which is balanced by a counter-propagating 
current carried by the unpaired quarks near their 
Fermi surfaces that are not in the pairing ribbons. 
Hence, the state carries no net current. 

After solving a gap equation for A and then eval- 
uating the free energy of the phase with condensate 
(|12|1 . one finds that this simplest "crystalline" phase 
is favored over two-flavor quark matter with either 
no pairing or BCS pairing only within a narrow win- 
dow 

0.707 A2SC <S^i< 0.754 Ajsc , (13) 

where A2SC is the gap parameter for the two-flavor 
phase with 2SC (2-flavor, 2-color) BCS pairing found 
at 5 fi ^ 0. At the upper boundary of this window, 
A — > and one finds a second order phase tran- 
sition between the crystalline and unpaired phases. 
At the lower boundary, there is a first order transi- 
tion between the crystalline and BCS paired phases. 
The crystalline phase persists in the weak coupling 
limit only if 6fi/A2sc is held fixed, within the win- 
dow H13() . while the standard weak-coupling limit 
A2Sc/m ^ is taken. Looking ahead to our con- 
text, and recalling that in three-flavor quark matter 
Sfi = Mg/(8/^), we see that at high densities one 
finds the CFL phase (which is the three-flavor quark 
matter BCS phase) and in some window of lower 
densities one finds a crystalline phase. In the vicinity 
of the second order transition, where A — ^ and in 
particular where A/Sfj, and, consequently given 
((T^ . A/A2SC — > a Ginzburg-Landau expansion of 
the free energy order by order in powers of A is con- 
trolled. Analysis within an NJL model shows that 
the results for A((5/x) become accurate in the limit 
6fi 0.754 A2SC where A — > 0, as must be the case, 
and show that the Ginzburg-Landau approximation 
underestimates A{Sfi) at all Sfi |3(]ll35l|. 

The Ginzburg-Landau analysis can then be ap- 
plied to more complicated crystal structures in which 
Cooper pairs with several different q's, all with 
the same l eng th but pointing in different direc- 
tions, arise p35l |. This analysis indicates that a face- 
centered cubic structure constructed as the sum of 



eight plane waves with q's pointing at the corners of 
a cube is favored, but it does not permit a quanti- 
tative evaluation of A((5/i). The Ginzburg-Landau 
expansion of the free energy has terms that are 
quartic and sextic in A whose coefficients are both 
large in magnitude and negative. To this order, J7 
is not bounded from below. This means that the 
Ginzburg-Landau analysis predicts a strong first or- 
der phase transition between the crystalline and un- 
paired phase, at some Sfi significantly larger than 
0.754 A2SC1 meaning that the crystalline phase oc- 
curs over a range of dfi that is much wider than 
but it precludes the quantitative evaluation of the Sfx 
at which the transition occurs, of A, or of fi. 

We shall find that in three-flavor quark mat- 
ter, all the crystalline phases that we analyze have 
Ginzburg-Landau free energies with positive sextic 
coefficient, meaning that they can be used to eval- 
uate A, n and the location of the transition from 
unpaired quark matter to the crystalline phase with 
a postulated crystal structure. For the most favored 
crystal structures, we find that the window in pa- 
rameter space in which they occur is given by 
which is in no sense narrow. 



D. Crystalline color superconductivity in 
neutral three-flavor quark matter 

Our purpose in this paper is to analyze three- 
flavor crystalline color superconductivity, with con- 
densates as in Eq. |^ for a variety of choices of the 
sets of Qj's and q^'s, i.e. for a variety of crystal 
structures. We shall make weak coupling (namely 
AojSfi ^ fi) and Ginzburg-Landau (namely A <^ 
AQ,Sfi) approximations througout. 

The analysis of neutrality in three-flavor quark 
matter in a crystalline color superconducting phase 
is very simple in the Ginzburg-Landau limit in which 
A <^ Sfx: because the construction of this phase does 
not involve rearranging any Fermi momenta prior to 
pairing, and because the assumption A ^ Sfi im- 
plies that the pairing does not significantly change 
any number densities, neutrality is achieved with 
the same chemical potentials iie = M'^ / (4/i) and 
fJ-s = Ms = as in unpaired quark matter, and with 
Fermi momenta given in Eqs. JB)), and (|10() as in 
unpaired quark matter. This result is correct only 
in the Ginzburg-Landau limit. 
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We consider a condensate of the form 



{lpia{.x)C^5'4)jl3{x)) (X 

1=1 q?e{qj 



J2 A/e2*'i° '-e/„^e/,, , 



(14) 

where q°, and q§ and Ai, A2 and A3 are the 
wave vectors and gap parameters describing pair- 
ing between the (d, s), {u,s) and {u,d) quarks re- 
spectively, whose Fermi momenta are split by 
26iJ,2 and 26^3 respectively. From (|1U|) . we see that 
Sfj,2 = <5/^3 = Sni/2 = Mf/{8fi). For each /, {q^} is 
a set of momentum vectors that define the periodic 
spatial modulation of the crystalline condensate de- 
scribing pairing between the quarks whose flavor is 
not /, and whose color is not /. Our goal in this pa- 
per is to compare condensates with different choices 
of {q^l's, that is with different crystal structures. 
To shorten expressions, we will henceforth write 



(15) 



qjefq/} 



The condensate (|14|l has the color-flavor structure of 
the CFL condensate (obtained by setting all q's to 
zero) and is the natural generalization to nontrivial 
crystal structures of the condensate previously ana- 
lyzed in Refs. |39l.l4ll|. in which each {qj} contained 
only a single vector. 

In the derivation of the Ginzburg-Landau approx- 
imation in Section IV, we shall make no further as- 
sumptions. However, in Sections V and VI when we 
evaluate the Ginzburg-Landau coefficients and give 
our results, we shall make the further simplifying 
assumption that Ai = 0. Given that S^i is twice 
(5/i2 or Sfi3, it seems reasonable that Ai ^ A2, A3. 
We leave a quantitative investigation of condensates 
with Ai 7^ to future work. 



E. NJL Model, and Mean-Field 
Approximation 

As discussed in Section I, we shall work in an NJL 
model in which the quarks interact via a point-like 
four-quark interaction, with the quantum numbers 
of single- gluon exchange, analyzed in mean field the- 
ory. By this we mean that the interaction term 
added to the Lagrangian Q is 



where we have suppressed the color and flavor in- 
dices that we showed explicitly in ^ , and have con- 
tinued to suppress the Dirac indices. The full ex- 
pression for r^" is (T^'')ai,i33 = 7'^iT^)a0Sij, where 
the are the color Gell-Mann matrices. The NJL 
coupling constant A has dimension -2, meaning that 
an ultraviolet cutoff A must be introduced as a sec- 
ond parameter in order to fully specify the interac- 
tion. Defining A as the restriction that momentum 
integrals be restricted to a shell around the Fermi 
surface, ^ — A < |p| <^-|-A, the CFL gap parame- 
ter can then be evaluated: 0, ISSj 



2SAcxp 



^2 1 



2/i2A 



(17) 



We shall see in subsequent sections that in the limit 
in which which A <^ Ao,Sfi <C /i, all our results 
can be expressed in terms of Aq; neither A nor A 
shall appear. This reflects the fact that in this limit 
the physics of interest is dominated by quarks near 
the Fermi surfaces, not near A, and so once Aq is 
used as the parameter describing the strength of the 
attraction between quarks, A is no longer visible; 
the cutoff A only appears in the relation between 
Aq and A, not in any comparison among different 
possible paired phases. In our numerical evaluations 
in Section VI, we shall take fj, = 500 MeV, A = 
100 MeV, and adjust A to be such that Ao is 25 MeV. 

In the mean-field approximation, the interaction 
Lagrangian (|16|l takes on the form 

Antcraction = ^ijA{x)i:^ + ^l/j'^ A{x)lP , (18) 

where A(a;) is related to the diquark condensate by 
the relations 



A{x) = ^Ar^-(^v^^)(r^,)^ 
A{x) = ^A(r^'')^(^^^)r^, 

= 7°At(a;)7° . 



(19) 



-'interaction 



;A(v;r^'^^)(^r^,^) 



(16) 



The ansatz (|14ll can now be made precise: we take 
A{x) = Acf{x) <E) C-f^ , (20) 

with 

3 

Acf(x)„,,^, =^^A(q?)e2*'J°-^e,„^e,,, . (21) 

1=1 q? 
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We have introduced notation that allows for the 
possibility of gap parameters A(qj) with different 
magnitudes for different / and for different a. In 
fact, we shall only consider circumstances in which 
A(qf ) = A/, as in H14|l . However, it will be very con- 
venient in subsequent sections to keep track of which 
A/ in a comphcated equation "goes with" which qj, 
making this notation useful. 

The full Lagrangian, given by the sum of (jj)) and 
(|18|l . is then quadratic and can be written very sim- 
ply upon introducing the two component Nambu- 
Gorkov spinor 



X 



and hence x — ( "0 "0"^ ) 



in terms of which 



l_ i^ + l/i A(x) 
2^ I A(x) {i^-iliY ' ^ 



(22) 



(23) 



Here, i/l = /i7o and ^, is the matrix Q), which we 
have argued simplifies to 



^ = (5"'^ (g) diag ^ld, Ms 



(24) 



with the flavor chemical potentials given simply by 
l|l(J|l . In subsequent sections, we shall also often use 
the notation j^i = /^i7o, with i = 1,2,3 correspond- 
ing to M, d, s respectively. 

The propagator corresponding to the Lagrangian 
is given by 

iG{x,x') iF{x,x') 
iF{x,x') iG(x,x') 

(25) 

where G and G are the "normal" components of the 
propagator and F and F are the "anomalous" com- 
ponents. They satisfy the coupled differential equa- 
tions 



( 



i^_ + t/L A(a-) \(G{x.x') F{x,x') 
A{x) {i^- i^f )[ F{x,x') G{x,x') 



= (J l)s^'Hx~x'). 

We can now rewrite (|19ll as 

Aix) = jXr^''Fix,x)iTA.f 

A{x) = ^\(r^''fF{x,x)rA. , 



(26) 



(27) 



either one of which is the self-consistency equation, 
or gap equation, that we must solve. Without fur- 
ther approximation, (|27|) is not tractable. It yields 
an infinite set of coupled gap equations, one for each 
A(qj), because without further approximation it is 
not consistent to choose finite sets {q/}. When sev- 
eral plane waves are present in the condensate, they 
induce an infinite tower of higher momentum con- 
densates |35j . The reason why the Ginzburg-Landau 
approximation, to which we now turn, is such a sim- 
plification is that it eliminates these higher harmon- 
ics. 



III. GINZBURG-LANDAU 
APPROXIMATION: INTRODUCTION 

The form of the Ginzburg-Landau expansion of 
the free energy can be derived using only general ar- 
guments. This, combined with results for two-flavor 
crystalline color superconductivity from Ref. |35| . 
will allow us to draw some partial conclusions in this 
Section. 

We shall only consider crystal structures in which 
all the vectors q^ in the crystal structure {q/} are 
"equivalent" . By this we mean that a rigid rotation 
of the crystal structure can be found which maps any 
q" to any other q^ leaving the set {q/} invariant. 
For such crystal structures, A(qj) = A/, meaning 
that the free energy is a function only of Ai, A2 and 
A3. As explained in Section II. D, the chemical po- 
tentials that maintain neutrality in three-fiavor crys- 
talline color superconducting quark matter are the 
same as those in neutral unpaired three-fiavor quark 
matter. Therefore, 



crystalline 



= n 



unpaired 



f7(Ai,A2,A3) , (28) 



with Ounpaired givcu in {Tj) with JHJ, and with 
ri(0,0, 0) = 0. Our task is to evaluate the conden- 
sation energy r2(Ai, A2, A3). Since our Lagrangian 
is baryon number conserving and contains no weak 
interactions, it is invariant under a global U{1) sym- 
metry for each flavor. This means that Q must be 
invariant under A/ e"^^A/ for each /, meaning 
that each of the three A/'s can only appear in the 
combination AJA/. (Of course, the ground state 
can and does break these U{1) symmetries sponta- 
neously; what we need in the argument we are mak- 
ing here is only that they are not explicitly broken 
in the Lagrangian.) We conclude that if we expand 
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51(Ai, A2, A3) in powers of the A/'s up to sextic or- der, it must take the form 



I 

\ {^Pi{A*Aif + |] /3„ AJA,A}A^j 

i 7/(AJ A,)3 + AJ A,A}AjA}Aj + 7123 AJ AiA^AaA^Ag j 



(29) 



where we have made various notational choices for 
later convenience. The overall prefactor of 2/i^/7r^ 
is the density of states at the Fermi surface of un- 
paired quark matter with Ms = 0; it will prove con- 
venient that we have defined all the coefficients in 
the Ginzburg-Landau expansion of the free energy 
relative to this. We have defined Pj = dim{q/}, 
the number of plane waves in the crystal struc- 
ture for the condensate describing pairing between 
quarks whose flavor and color are not /. Writ- 
ing the prefactor P/ multiplying the quadratic term 
and writing the factors of ^ and multiplying the 
quartic and sextic terms ensures that the a/, /?/ 
and c oefficients are defined the same way as in 
Ref. [33. The form of the Ginzburg-Landau ex- 
pansion (|29|l is model independent, whereas the ex- 
pressions for the coefficients a/, /?/, /3/j, 77, 7/jj 
and 7123 for a given ansatz for the crystal structure 
are model-dependent. In Section IV we shall derive 
the Ginzburg-Landau approximation to our model, 
yielding expressions for these coefiicients which we 
then evaluate in Section V. 



We see in Eq. H29() that there are some coefficients 
— namely a/, /?/ and 7/ — which multiply polyno- 
mials involving only a single A/. Suppose that wc 
keep a single A/ nonzero, setting the other two to 
zero. This reduces the problem to one with two- 
flavor pairing only, and the Ginzburg-Landau co- 
efficients for this problem have been calculated for 
many different crystal structures in Ref. [s^. We 
can then immediately use these coefficients, called 
a, (3 and 7 in Ref. |35| . to determine our a/, /3/ and 



7/. Using ai as an example, wc conclude that 



ai = a{qi,6fii) 



2qi 



qi + Sfii 
qi - 5iii 



A2 

^2SC 



4(g| - 5^i]) 



(30) 



where is the splitting between the Fermi surfaces 
of the quarks with the two flavors other than / and 
qi = |qj| is the length of the q- vectors in the set 
{q/}. (We shall see momentarily why all have the 
same length.) In (|^ . A2SC is the gap parameter 
in the BCS state obtained with S^ii = and A/ 
nonzero with the other two gap parameters set to 
zero. Assuming that Aq <C ^, this gap parameter for 
2SC (2-flavor, 2-color) BCS pairing is given by 0,113 



A2SC = 23Ao 



(31) 



For a given S^ij and Ao, aj given in (|30|l is minimized 
when 111, 35, 42J 



qi ~ rjSj-ii with ry = 1.1997 . 



(32) 



In the Ginzburg-Landau approximation, in which 
the A/ are assumed to be small, we must first min- 
imize the quadratic contribution to the free energy, 
before proceeding to investigate the consequences 
of the quartic and sextic contributions. Minimiz- 
ing aj fixes the length of all the q-vectors in the 
set {q/}, thus eliminating the possibility of higher 
harmonics. It is helpful to imagine the (three) sets 
{qj} as representing the vertices of (three) polyhe- 
dra in momentum space. By minimizing a/, we have 
learned that each polyhedron {q^} can be inscribed 
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in a sphere of radius "qSfii. From the quadratic 
contribution to the free energy, we do not learn 
anything about what shape polyhedra arc prefer- 
able. In fact, the quadratic analysis in isolation 
would indicate that if a/ < (which happens for 
S^i < 0.754 A2sc) then modes with arbitarily many 
different q/'s should condense. It is the quartic 
and sextic coefficients that describe the interaction 
among the modes, and hence control what shape 
polyhedra are in fact preferable. 

The quartic and sextic coefficients /?/ and 7/ can 
also be taken directly from the two-flavor results of 
Ref. • They are given by P/S^'^ and "f/S^j where 
j3 and 7 are dimensionless quantities depending only 
on the directions of the vectors in the set {q^}. They 
have been evaluated for many crystal structures in 
Ref. [35I I. resulting in two qualitative conclusions. 
Recall that, as reviewed in Section II. C, the presence 
of a condensate with some qf corresponds to pairing 
on a ring on each Fermi surface with opening angle 
67.1°. The first qualitative conclusion is that any 
crystal structure in which there are two q^'s whose 
pairing rings intersect has very large, positive, val- 
ues of both Pj and 7/, meaning that it is strongly 
disfavored. The second conclusion is that regular 
structures, those in which there are many ways of 
adding four or six qj's to form closed figures in mo- 
mentum space, are favored. Consequently, according 
to Ref. j35l | the favored crystal structure in the two- 
flavor case has 8 qj's pointing towards the corners 
of a cube. Choosing the polyhedron in momentum 
space to be a cube yields a face-centered cubic mod- 
ulation of the condensate in position space. 

Because the /?/ and 7/ coefficients in our problem 
can be taken over directly from the two-flavor anal- 
ysis, we can expect that it will be unfavorable for 
any of the three sets {qj} to have more than eight 
vectors, or to have any vectors closer together than 
67.1°. At this point we cannot exclude the possi- 
bility that the large positive /3/ and 7/ indicating 
an unfavorable {q/} could be offset by large nega- 
tive values for the other coefficients which we cannot 
read off from the two-flavor analysis. However, what 
we shall instead find in Section VI is that /3/j and 
7//J are positive in all cases that we have investi- 
gated. This means that we know of no exceptions to 
the rule that if a particular {q^} is unfavorable as 
a two-flavor crystal structure, then any three-flavor 
condensate in which this set of q-vectors describes 
either the Ai, A2 or A3 crystal structure is also dis- 
favored. 



In Section IV we shall use our microscopic model 
to derive expressions for all the coefficients in the 
Ginzburg-Landau expansion H29|) , including rederiv- 
ing those which we have taken above from the two- 
flavor analysis of Ref. [s^ . The coefficients that we 
cannot simply read off from a two-flavor analysis are 
those that multiply terms involving more than one 
A/ and hence describe the interaction between the 
three different A/'s. Before evaluating the expres- 
sions for the coefficients in Section V, we shall make 
the further simpifying assumption that Ai = 0, be- 
cause the separation Sfii between the d and s Fermi 
surfaces is twice as large as that between either and 
the intervening u Fermi surface. This simplifies (I29|l 
considerably, eliminating the 7123 term and all the 
Pij and 7/7 J terms except ^32, 7223 and 7332. 

IV. THE GINZBURG-LANDAU 
APPROXIMATION: DERIVATION 

We now derive the Ginburg-Landau approxima- 
tion to the NJL model specified in Section II. We 
proceed by first making a Ginzburg-Landau approx- 
imation to the gap equation, and then formally in- 
tegrate the gap equation in order to obtain the free 
energy, since the gap equation is the variation of the 
free energy with respect to the gap parameters. 

The gap equation (|27|l with which we closed Sec- 
tion II is an infinite set of coupled equations, one 
for each A(qj), with each equation containing arbi- 
trarily high powers of the A's. In order to make a 
Ginzburg-Landau expansion, order by order in pow- 
ers of the A's, we first integrate (|26|l . obtaining 

G(„0^C<»),.,.V/^'.G'»'(.,.)A(,)F(=,.') 

F{x, 2;') = - y" G'"' {x, z)A{z)G{z, x') 

(33) 

with = [i^ + if)-^ and GW = {{i^ - 

We then expand these equations order by order in 

A(.t) by iterating them. To fifth order, for F wc find 

F = - g(°)ag(") - g(")ag(")agWag(°) 

_ (5(0)^^(0)^^(0)^^(0) ^^(0)^^(0) ^0(^7) ^ 

(34) 

where we have suppressed space-time coordinates 
and integrals for simplicity. Wc then substitute this 
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7At 



FIG. 1: The gap equation. The labels a, (3 represent the external color indices and i, j represent the external flavor 
indices. All the other color-flavor indices are contracted. Acf (and A^^) are matrices of the form 1211 and carry the 
same color and flavor indices as the neighbouring propagators. The dashed lines represent the propagator {i^ — 
and the solid lines represent {i<^ + i/l)^^ . Evaluating the gap equation involves substituting II2H for Acf, doing the 
contraction over the internal color-flavor indices, and evaluating the loop integrals in momentum space. 



expansion for F into the right-hand side of the gap 
equation for A(a;) in lf77|) . After using the C75 
Dirac structure of our ansatz (|20|l and the identity 



-7^ to simplify the expression, we 



obtain the gap equation satisfied by 1^cf{x)^ the 
part of our ansatz H2UI21|) that describes the color, 
flavor and spatial form of our condensate. To order 
A^, we find 



7A 



1 



tA 



CF 



CF , 



tA 



CF 



-A^ — 



1 



(35) 



7A 



CF 



tA 



CF 



1 t 1 ^ 1 1 



{x.x) 



where the differential operators act on everything to 
their right and where we have continued to simplify 
the notation by not writing the space-time, color and 
flavor arguments of the Acf^s and by not writing the 
integrals. We then use the color Fierz identity 

(i")a'a(i")/3'/3 - {-^S^^,dpp> + 26a,p6a.p') (36) 

to rewrite (|35|l as shown diagrammatically in Fig. ^ 
As written in H35|) and shown in Fig. what oc- 
curs on the left-hand side of the gap equation is the 



space-dependent condensate from (|21() . 

AJ,^(x)„,^, =^6,„^6,,,^A*(q?)e-2»'i?- , 

(37) 

whereas we now wish to turn the gap equation into 
a set of coupled equations for the constants A(qJ). 
Doing so requires simplification of the color-flavor 
structure of the right-hand side. Our ansatz for the 
color-flavor structure of the condensate, on the left- 
hand side, is antisymmetric in both color and flavor. 
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However, direct evaluation of the right-hand side 
yields terms that are symmetric in color and flavor, 
in addition to the desired terms that are antisym- 
metric in both. This circumstance is familiar from 
the analysis of the CFL phase P, , whose color- 
flavor structure we are after all employing. In the 
presence of a color and flavor antisymmetric conden- 
sate, a symmetric condensate must also be generated 
because doing so does not change any symmetries. 
The same argument applies here also. In the CFL 
phase, the symmetric condensate is both quantita- 
tively and paramctrically suppressed relative to the 
antisymmetric condensate, which is understandable 
based on the basic fact that the QCD interaction is 
attractive in the antisymmetric channel and repul- 
sive in the symmetric channel. We therefore expect 
that here too if we were to include color and flavor 
symmetric condensates in our ansatz and solve for 
them, they would prove to be suppressed relative to 
the antisymmetric condensates, and furthermore ex- 
pect that, as in the CFL phase, their inclusion would 
have negligible impact on the value of the dominant 



antisymmetric condensate. Hence, we drop the color 
and flavor symmetric terms occurring on the right- 
hand side of the gap equation. Upon so doing, the 
right-hand side of the gap equation, which we denote 
Rai,i3j, has the structure 

Raz,i3j{x) = ^ Ri{r)eiapeiij (38) 

Because eiap^Hj are linearly independent tensors for 
each value of /, in order for the gap equation to be 
satisfied for all values of a, /3, i and j we must have 

^A*(q?)e-2''^^"■ = i?z(r) (39) 
q? 

for all three values of /. This is a set of Yli Pi 
coupled equations for the undetermined constants 
A(qj). (Recall that Pi is the number of vectors 
in the set {q^}.) After transforming to momentum 
space, these gap equations can be written as follows: 
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+ ^ A*(qf)A(qf)A*(qf)J,,,,(qJ,qJ,qf,q?)5q5„q.+q._q„ 
+ E A*(qj)A(q})A*(qf)J)c/fc,/(q^„q},q?,q?)^q._qc+q._q. 

./ ql^qfjq? 

q?qfq?qjqf 

^ A*(q^^)A(q})A*(q^)A(ql)A*(q{)Wj.j(q^^,q},q?,q^q{,q?)'5<^._q.+q._^ 
qjqjq^q/q^ 

4E E A*(q'j)A(q})A*(q^,)A(q})A*(qf)W,.,(q^,,q},q^,,q},q{,q?)J^._^.^^._^^^ 
q'qjqjqjq^ 



,/_qa 



c / 



-i 5] ^ A*(q^,)A(q})A*(qf )A(q^^)A*(q^) /C™//,/(q'7, q}, q?, q?,, q^, q?) ^^.^c 
-\Y. H A*(q^^)A(q^)A*(q^)A(q})A*(q^)/CK/Ji^7j(q^j,q^,q?,q},q^,q?)<5q._q^+q._q. 



-q}+qi<-q? 



(40) 



where we have introduced a lot of notation that we 
now define and explain. First, recall from H15() that 
means a sum over all the q^'s in the set {q^}. 
The (5's are therefore Kronecker (5's, indicating that 
only those q- vectors that can be arranged to form a 
certain closed two-, four- or six-sided figure in mo- 
mentum space are to be included in the sum. The 
sums over J are always understood to be sums over 
J ^ I, and the sums over K are always understood 
to be sums over K ^ J and K ^ I . The remain- 
ing flavor subscripts in some terms which are not 



summed, denoted j or fc, must always be chosen not 
equal to each other, not equal to /, and not equal to 
J if J occurs. (This appears to leave an ambiguity 
related to the exchange of j and k in terms where 
both occur, but we shall see that the functions IT, 
S and K, each have a cyclic symmetry that ensures 
that the two apparent choices of j and k are equiva- 
lent.) The functions H, J and JC are proportional to 
the various loop integrals that appear in the evalu- 
ation of the Feynman diagrams in the gap equation 
of Fig. n] They are given by 
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nij(ki,k2) = — 



7m 



Jij,fc,i(ki,k2,k3,k4) 



1 

7m 



+ 2}ti - 2% + 2^^: 



■3 m) 



A;ij,fc,;,m,n(ki,k2,k3,k4,k5,k6) = 5-7^ 



1 



^2 ' J {2^Y{i,- i^,){i> + 2]^i+ilL,){i) + 2]/ii-2yi2 + i^k) 

1 

1 

r^7M ' 



(j) + 2]^i - 2]/i2 + 2^^3 - 21^4 + 2^^5 



(41) 



where i/li = 7°/Xi and = (0, k)^j7^ = — k • 7. The 
subscripts i, j etc. on the functions 11, J and /C 
are flavor indices that give the flavor of the quark 
hnes in the propagators going around the loops in 
Fig. n In each term in H4U|I the choice of flavor in- 
dices in n, or /C is determined by the requirement 
that a given A(qj) must connect two propagators 
for quarks with flavors different from each other and 
/. For example, A3 always connects a u and a d 
quark. The easiest way to see how this provides the 
explanation for the (perhaps initially peculiar look- 
ing) prescriptions for the and /C functions in each 
term in the gap equations H4(J|) is to examine Fig. |21 
below, which depicts examples of the contributions 



of n, J and K. to the free energy which we shall 
discuss next. 

The gap equations that we have derived must be 
equivalent to the set of equations 9ri/9A(q°) = 0, 
because solutions to the gap equation are station- 
ary points of the free energy fi. This means that 
integrating the gap equations determines 17 up to an 
overall multiplicative constant, which we can fix by 
requiring that we reproduce known results for the 
single-plane wave condensates, and up to an addi- 
tive constant which we fix by the requirement that 

^^crystalline = f^unpaircd whcu all A(qJ) are SCt tO 

zero. We find 
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^^({A(q?)}) = ^E 



^ A*(qJ)A(q?)('n,,(q?,qJ) + -^' 



2 / '^q?-q? 



+ 1 A*(qJ)A(qJ)A*(qf)A(q?)^,fc,fe(qJ,qf,qf,q?)5q._qe+q._q. 
+ ^ E E A*(q!})A(q})A*(q?)A(q?) Jkikjici'j, q}, q?, q?)(5q. 

.J>I q'',q;jq?q? 

i ^ A*(qJ)A(qJ)A*(q?)A(q|)A*(qf)A(q?)/C,,,,,,(q^q5,q^,ql,qf,q?)^^._^ 
qjq/q/qfq/q? 

i ^ 5] A*(q^^)A(q})A*(q?)A(qf)A*(q{)A(q?) 



'^'^-^qjqjq^q/q/q? 



ICkIkJkj{q^J, q}, q?, q?, > q/)'^q!j-q- +qd-q=+qf-c 



-q? 

1 5^ 5^ A*(q5)A(q})A*(q^)A(q^)A*(q5)A(q?) 

J5^A:#/7^./q5qCq^.q^q/qa 

/Cjif/,//A'(q/, q}, qi, qx, q^, q?)'5q5-qc +q^^,-qj^+q/-q? 
1 ^ ^ A*(q^^)A(q^)A*(qf)A(q})A*(q^)A(q?) 
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ICKiJKuiq!}, q^, qf , q}, q^, q?)'^q|}-q^+q?-q5+q{^-q? 



(42) 



As in (|40|) . in each term the flavor indices j and k (or just fc) that are not summed over are understood to 
differ from each other and from the summed indices / (or / and J). 

As wc discussed in Section III, we shall only consider crystal structures in which each of the three sets 
{q^^} are regular, in the sense that all the q° in one set {q^} are equivalent. This means that A(qj) = A/, 
which simplifies the free energy (|42|) to the form H29|l which we derived on general grounds in Section III 
and which we reproduce here 



0(Ai,A2,A3)=^ ^P,a/AJA, + i |^^/3,(AlA,)2 + ^^^/3,jA*A,A}A^j 

/ ' \n (43) 

+ \ ^i{A}Aif + 7/JjAJ A,A}AjA}Aj + 7i23AJ AiA* A2A5A3 
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^p + 2(q^-q§) g 



p + 2(q§-q§ + qg) 




P + 2(q^-q§) 



A2K P + 2(q§-q§ + q^) 



p + 2(q^-q§) ^,^<§)-^p + 2(q^-q§ + q^) 



P + 2q§ 




P + 2(q^ - q§) , , -(A|)-^j; + 2(q^ - q§ + qf ) 



p + 2(q^-q§) 



/+2(q^-q§) 



+2(q§-q§ + q^) 



P + 2q^ 




p + 2(q§ - q§ 
2(qf-qi) 



- /p + 2(q^-q§) 

d „6 _L J\ 



+2(q^-q^ + q2M 



FIG. 2: Examples of contributions to the free energy. The five diagrams depict a Tlud contribution to asAsAs, a 
Judv.d contribution to /33(A|A3)^, a J^dua contribution to /332 A3 A3 A2 A2, a ICududud contribution to 73(A3A3)^ and 
a ICudusua contribution to 7322A|A3(A2A2)^. 

for continuity. Now. however, we have obtained expUcit expressions for all of the coefficients: 
a, =n,fe(q?,q?) 



J,-;.,A:(q^q/,q?,q?)^q;- 



P. 



JI 



7123 



E J'fc/fc./(q'j,q},q?,q?)'5q.-q5- 
qjq?q? 

5Z iCjkokjk (q/, q/, qf , qf , q{q?)'5q5_q^+q^ 



-qj+qi--q/ 



q^q^q^q^q^q? 



(44) 



X! ^fc/fc jfc J (qj , q} , q/ , q? , q{ , q?)'^q'', -q^+q^ 



qjqf/qjqjqfq? 



qf+qj-q} 



XI ^jif/j/K(q/,qj,qi-,qif,qj,q/)'5q6_ 



q;-q5+q?f-qA-+q7-q? 



7 E E ^™^/j(qj:qK,q?,q},qt,q?)'5q._qc^,+q._qe+q^_qo . 



I^J^K^I qt qc^^qdqe q/ q. 
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Here again, the unsummcd indices j and k are cho- 
sen as described previously. Since the free energy 
(|48|l is invariant under phase rotations of the A/ we 
can henceforth take all the A/ real and positive. In 
Fig. 13 we give examples of contributions to the free 
energy. These examples should make clear the choice 
of flavor subscripts on the J''s and /C's in H44|) and 
consequently in H4()|l . They also illustrate the origin 
of the Kronecker 5's in so many of the expressions in 
this section: each insertion of a A(qj) (or A*(qf)) 
adds (or subtracts) momentum 2qf to (from) the 
loop, meaning that the Kronecker (5's arise due to 
momentum conservation. The diagrams also illus- 
trate that n, and /C are invariant under simulta- 
neous cyclic permutation of their flavor indices and 
momentum arguments, as this corresponds simply 
to rotating the corresponding diagrams. 

We have succeeded in deriving expressions for the 
Ginzburg-Landau coefficients in our model; we shall 
turn to evaluating them in the next section. Re- 
call, however, that upon setting Ai = and keep- 
ing in mind that we can obtain results for or, /3/ 
and 7/ from the two-flavor analyses in Ref. [23 j 
that we need to do is evaluate /332, 7233 and 7322 
for the crystal structures we wish to investigate. We 
shall largely focus on crystal structures for which 
{€[2} and {q3} are "exchange symmetric", meaning 
that there is a sequence of rigid rotations and reflec- 
tions which when applied to all the vectors in {q2} 
and {q3} together has the effect of exchanging {(I2} 
and {q3}. If we choose an exchange symmetric crys- 
tal structure, upon making the approximation that 
(5/i2 = and restricting our attention to solutions 
with A2 = A3 we have the further simpliflcation 
that 7322 = 7233- Once we learn how to evaluate 
the loop integrals and JC in the next Section, we 
will then in Section VI evaluate P32 and 7322 for var- 
ious crystal structures, enabling us to evaluate the 
magnitudes of their gaps and condensation energies. 



ming JTudus (Qs, qi, q2i I2) O'^^r all those vectors q| 
and q§ in the set {q3} and all those vectors q2 and 
q2 in the set {q2} which satisfy q^ — q3+q2 — q2 =0, 
forming a closed four-sided flgure in momentum 
space. Understanding how to evaluate the loop inte- 
grals n, J and K. requires some explanation, which is 
our goal in this section. Performing the sum required 
to evaluate a given Ginzburg-Landau coefficient is 
then just bookkeeping, albeit nontrivial bookkeep- 
ing for complicated crystal structures. 

We are working in a weak-coupling limit in which 
Sfi, |q| ~ q ^ V^l^i and A2SC are all much smaller 
than /i. This means that we can choose our cutoff 
A such that (5/^, g, A2SC ^ A ^ Because A <^ fj,, 
the integration measure in the expressions (|41|l for 
n, J' and /C simplifies as follows: 



-A 2 



dp 

4^ ' 



(45) 



where s = \p\ — fi. We now see by power counting 
that n is log-divergent as we take A ^ Sfi, q, A2SC 
whereas both and /C are A-independent in the 
large A limit. Thus, in evaluating J7 and /C, we can 
safely take A ^ 00 whereas we must keep A in the 
problem for a little longer in analyzing 11. Explicit 
evaluation of 11 yields 



n„d(q3,q3) 



1 + log J— 

293 \q3-(>fj,3 



(46) 



We can now use 



A. 



V. CALCULATING GINZBURG-LANDAU 
COEFFICIENTS 

Calculating the Ginzburg-Landau coefficients (|44|l 
that specify ri(Ai, A2, A3) for a given crystal struc- 
ture involves first evaluating the loop integrals 11, 

and /C, defined in H41|) . and then summing those 
that contribute to a given Ginzburg-Landau coef- 
ficient. For example, we see from (O that the 
Ginzburg-Landau coefficient /332 is given by sum- 



2SC 



2Ae 217^ 



(47) 



and the relation between and 11^^ given in 
to evaluate 03, obtaining the result (|30|l . No- 
tice that a/ depends on A and A only through 
A2SC! and depends only on the ratios <7//A2sc and 
S fii / A2SC- As discussed in Section III, a/ is nega- 
tive for Sfii /A2SC < 0.754, and for a given value of 
this ratio for which a/ < 0, a/ is most negative for 
g7/A2sc = '7'5m//A2sc with 77 — 1.1997. We there- 
fore set qi = rj S^i henceforth and upon so doing 
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obtain 

where in the last hne we have used the definition of 
rj derived from H30() . 

The evaluation of f3i and 7/ is described in 
Ref. • From the integration measure (|45|l and 
the definitions of J' and K. (|41|l we see that /3/ and 
7/ have dimension -2 and -4, respectively. Since they 
are independent of A as long as A Sfi, q, A2SC, and 
since A nowhere appears in their definition, there is 
no need to introduce A2SC • This means that the only 
dimensionful quantity on which they can depend is 
6fii (since qi = rjSiij and since the propagators are 
independent of fj. in the weak-coupling limit) and so 
we can write 

/?, = |L and 11 = j^, (49) 
I 



where Pi and 77 are dimensionless quantities that 
depend only on the shape of the polyhedron de- 
scribed by the set of vectors {q/}. The evaluation 
of the J and /C loop integrals occurring in (3 and 7 
is described in Ref. [s^, and results for many two- 
flavor crystal structures {q^} are tabulated there. 
The evaluation is similar to but simpler than the 
evaluation of (332 and 7322 , to which we now turn. 

P32 is the sum of Judus{(ll,q^,qi,q.2), where the 
momentum vectors satisfy 



qa - qs + q2 - q2 - . (5o) 



We now utilize the fact that |q|| = |q§| = 77(5 /is and 
|q2l = |q2l = ??i5/i2 where 6^3 and (5/i2 are similar 
in magnitude, but not precisely equal. (Recall from 
Section II. A that both are given by M^/(8/z) to this 
order, but that they differ at order Mf/ fj?.) Because 
5 7^ ^A'S: the condition (|^ can only be satisfied if 
qls = q-ji and q2 = q2. We must therefore evaluate 



J«dn.(q^,q^,q5,q^) = -^7^ 



(2^)4 _ ^.,){i> + 2^^ + ,/i,i){i> - ,/iu){i> + + ^) 



7m 



(51) 



We now expand the propagators in the weak-coupling limit, in which p°, s, |q|, (/i^ — //„) and (//„ — ^s) are 
all small compared to as follows: 

1 _ (pO + /z07° - (p + 2q) • 7 



(pO + /i» - |P + 2q|)(pO + M« + |P + 2q|) 

Mu7° - P ■ 7 

{p° + /i« - (m« - A^i) - IpI - 2q • p)(2^„) 



7" - p • 7 



(52) 



2 Vp° - s + - fJ.u) - 2q • p 



Similarly, 



7 + P • 7 



pO + s - (^i - fiu) + 2q • p 



(53) 



Eq. H51() then simplifies to 

J«...(q3,q3,q2,q2)-y ^ ^ ds 



(pO + s)2(p0 - s - p • 2q!; + 2dfi3){p" - s - p • 2q^ - 2Sfi2 



r 



where we have used (5 /is = 5 (/id — /in) and (5/i2 = ^(/"m ^ Ms)- 
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To integrate H54|) . we Wick rotate to ip* and 
then do the s integral by contour integration. This 
gives two contibutions with different sign factors, 
sign(p'*) and sign(— p''), which are complex conju- 
gates of each other. Combining the two, the inte- 
gration ovcrp* is of form 2^e /^^°° dp^{...) where we 



have started the integration from the infinitesimal 
positive number e instead of zero, thus dcfinining the 
principal value of the integral. The integration over 
can now be carried out safely to obtain 



1 f dp 



4 ./ Att 
1 



3?e 



ie - p • q| + Sfi3){ie - p • qg - 5^2 



dp 

4djj,2Sn3 J 47r 



5Re 



(55) 



(ie - i]p ■ q§ + l){ie - ?/p • qf - 1)_ 



where ri = = From rotational symmetry 

it follows that the value of H55|) depends only on 
the angle between the momentum vectors q| and 
q2, which we denote by (j). We therefore define the 
dimensionlcss quantities 



^32(0) = Sfi2Sn3Judus{cil,q^,q2,^2) 



(56) 



and, correspondingly, 

P32 = (5/i2(5^3/332 



(57) 



J32 can be evaluated analytically by using Feynman 
parameters to simplify the integrand in H55() . The 
result is 



J32(0) = 



^VC0sW2) [y^^2 3i,,2(^/2)-l 



arctan 



V sin2(0/2) - 1 
r]cos{(f)/2) 



(58) 



This completes the evaluation of the loop integral JT" 
needed to calculate /332 for any crystal structure. We 
summarize the calculation by noting that for a given 
crystal structure. P32 depends only on the shape of 
the polyhcdra defined by {q2} and {q^} and on their 
relative orientation, depends on the Fermi surface 
splittings and S^2; and is obtained using (|57|l 
with 



P32 = '^32(Zq^q^) , (59) 



where Jz2{(j^) is given by 

We turn now to the evaluation of 7322 . From H44(l , 



7322 =2 X! ^ndusus(q3,q3>q2>q2>q2>q2) : 

c,|,q§, 

(60) 

and we again use the fact that the q3's and q2's do 
not have precisely the same length to conclude that 
the momentum vectors must satisfy both 

q' = qs (61) 

and 

q2 - q2 + q^ - q2 = . (62) 

In the following expressions, it is always understood 
that H62|l is satisfied although we will not complicate 
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equations by eliminating one of the qj's in favor of 
the other three. We can see without calculation that, 
unlike J , JC will not reduce to depending only on a 
single angle between two momentum vectors. It will 
depend on the shape made by the four vectors 
satisfying l|62|l . which can in fact be specified by two 



angles, as well as on the angles that specify the di- 
rection of q| relative to the shape made by the four 

q2's- 

The expression for JC is given in H41|) and can also 
be read off from the bottom right Feynman diagram 
in Fig. 13 It is given by 



Y 



1 



After simplifying the propagators using (|52|1 . we can rewrite equation H63|l as 

1 



(63) 



(pO + s)2(pO + s + p-2(q^-q^)) 



1 



(pO - s - p • 2q^ + 2dfi3)ip° - s - p • 2q^ - 2Sfi2){p" ~s-p- 2(q^ - q^ + q^) - 2dfi2) 



(64) 



Unlike in the evaluation of Cfudus , we are not able to do the s and p'^ integrals analytically without introducing 
Feynman parameters to simplify the integrand at this stage, before doing any of the integrals. We introduce 
one set of Feynman parameters, xi, X2, to collect denominators of the form p^ + s + .. and another set, 
2/1) y2, Us, to collect the denominators of form — s + ... This reduces the integral to 

„i 2 / 2 \ 1 3 / 3 \ 

/C„d„sus(q3,q3,q2,q2)q2)q2) = / YldXnS il-^Xn\ '[[dy„iS il - "^ymj 

•^0 n=l \ 71=1 I ''^ m=l \ rn=l / 



dp 

47r 



oo 7 r-\-oo 



dp' 

2TTi 



ds 



4(1 -X2) 



(pO + s + 2a:2P- [q^-qi])3 
1 



(65) 



n 3 



- s - 2p • [yiq| + 2/2q2 + y3(q2 - qi + qD] + yi2(5/x3 - 2/22(S/i2 - ?/32(S/i2 

We now perform the p^ and s integrations in (|65|l . following steps analogous to the integration arising in 
the expression for Judus- i-C Wick rotate p^ to ip'*, do the s integral by contour integration, add the two 
complex conjugate contributions thus obtained to write the integration over as 23?e dp'^{...) and then 
perform the integration over p'^. This gives us 

/C„d„s„s(q3, qg, q2, q^, <4^^2) = / f[dx,i ( 1 - X! ^" ) / 11 ^^/m ^ {'^ ^ ) 

•^0 n=l \ n=l / •'^ m=l \ m=l J 



dp 

An 



1- X2 



p • [a;2(q2 - q2) - ivi^t + V2fii + 2;3(q2 - q2 + q2))] + yi'^Ms - V26^J.2 - v-,5^l2 



(66) 
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Finally, we do the dp integral and obtain 



•^0 n=l \ n=l I "^0 \ 

|2:2(q2 - qi) - (yiqs + y2q2 + y3(q2 - qi + q2))P + 3[2/i(S/i3 - yi&i>^2 - vz^^iif 



(1-^2) 



(67) 



|a;2(q2 " qi) - (yiqs + y2q2 + 2/3(q2 - qi + q2))P - W^^^z - y2^M2 - y?.^[^i + ie] 



Noting that we can replace (\2 by ^-iid q3 by T?q3, we conclude that, as expected, ICudusus depends only 
upon the shape of the polyhedra defined by {q2} and {q3} and on the Fermi surface splittings 6113 and Sn2- 
We cannot simplify H67|) further for general 6^2, 5^^. However, if we now set 8^2 = ^M3 = ^M, which is 
corrected only at order M*//i^, we can then factor out the dependence on the Fermi surface splitting, since 
the only dimensionful quantity in the integrand is then (5/i. Defining, for 5/i2 = (5/^3 = 



/C„d«s„^(q3,q3:q2,q2,q2:q2) = ^-f^'322(q3,q3,q2,q2,q2,q2) , 

and using |q^| = TjSjii, for all the momentum vectors, we find that ^^322 is given by 



(68) 



1 



■fi^322(q3,q3,q2,q2>q2>q2) = i / rfa;2(i-a;2) 

where 



dyi 



dy2^e- vM^2,yi,y2)\' + S{l - 2y,f 



[r;2|a(x2,yi,y2)P-(l-22;i)2 



& = X2 



(q' - 42) - (yiqs + y2q^ + (l - yi - y2)(q^ - q^ + 4) 



(69) 
(70) 



For general arguments we were not able to do 
the integrals that remain in 169|) analytically 
and therefore evaluated it numerically. Since 
■^322(q|,q3,q2,q2;q2'q2) is the limit of the tunc- 
tion7?322(q|,qtq^,qi,q{,qi,e)ase^0,wenu- 
merically evaluated the integral appearing in (|69|l 
at four values of e and extrapolated (using a cubic 
polynomial to fit the values) to e = 0. Finally 

7322 = 7322(5^^'^ (71) 

is found by summing /•C322 evaluated with 
all possible choices of momentum vectors 

(q3,qiq2>qi>q2>q2) satisfying ^ and mul- 
tiplying this sum by 3/2. 

VI. RESULTS 

A. Generalities 

We shall assume that Ai = throughout this sec- 
tion. As described previously, this simplification is 



motivated by the fact that Ai describes the pairing 
of d and s quarks, whose Fermi surfaces are twice 
as far apart from each other as either is from that 
of the u quarks. We shall focus most of our at- 
tention on exchange symmetric crystal structures, 
as defined at the end of Section IV, in which the 
polyhedra defined by {q2} and {q3} are related 
by a rigid rotation. In Section VI. D we will dis- 
cuss one example in which {q2} and {q3} are not 
exchange symmetric, and we have evaluated oth- 
ers. However, as none that we have investigated 
prove to be favorable, we shall make the notational 
simplifications that come with assuming that {q2} 
and are exchange symmetric, as this implies 

a2 = a3 = a, P2 = dim{q2} = P3 = dimlqg} = P, 
(32 — 03 = P and 7322 = 7233- The final simplifica- 
tion we employ is to make the approximation that 
S1J2 = 8fj,3 = Sfj. = M^/{8fi). As described in Sec- 
tion II. A, this approximation is corrected by terms of 
order M^/fi^. Upon making all these simplifying as- 
sumptions and approximations, the free energy ()43|l 
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reduces to 

+ l^ili^l + A3) + 7322(A^A^ + A^A^)) 

(72) 

where /?, 7, P32 and 7322 are the dimensionless con- 
stants that we must calculate for each crystal struc- 
ture as described in Section V, and where the 6fi- 
dependence of a is given by Eq. I^Hl- 

In order to find the extrema of ri(A2,A3) in 
(A2, A3)-space, it is convenient to write (A2,A3) 
as -\/2(Aj. COS0, Ar sin0) in terms of which the free 
energy ifT^ is given by 



n{Ar,d) 



2PaiSti)Al + -^PK 



5p? 



V2(5Ai^ 



(^32 - 2/3) + 



2A^ 

^(7322 -37)) sin' 29 



(73) 



Because sin' (20) has extrema only at = 7r/4 and 
9 = 0, 7r/2, we see that extrema of r2(A2,A3) ci- 
ther have A2 = A3 = A, or have one of A2 and 
A3 vanishing. The latter class of extrema are two- 
flavor crystalline phases. We are interested in the 
solutions with A2 = A3 = A. The stability of these 
solutions relative to those with only one of the A's 
nonzero appears to be controlled by the sign of the 
factor that multiplies sin' 26* in (|73|l . However, we 
shall show in Appendix^that the three-flavor crys- 
talline phases that we construct, with A2 = A3 = A, 
are electrically neutral whereas the two-flavor solu- 
tions in which only one of the A's is nonzero are not. 
Setting A2 = A3 = A, the free energy becomes 



n{A) 



2n^ 



2Pa{5ii)A^ 



2511 



A6 



-7off 



where we have defined 



/3eff ^2p + /?32 

7off = 27-4- 27-322 



(74) 



(75) 



We have arrived at a familiar-looking sextic order 
Ginzburg-Landau free energy function, whose coef- 
ficients we will evaluate for specific crystal structures 
in VLB and VI. D. First, however, we review the 
physics described by this free energy depending on 
whether /3eff and 7eff are positive or negative. 

If /?off and 7eff are both positive, the free energy 
l(7^ describes a second order phase transition be- 
tween the crystalline color superconducting phase 
and the normal phase at the Sfi at which a{Sii) 
changes sign. From (|48|l . this critical point occurs 
where (5/i = 0.754 A2sc- In plotting our results, we 
will take the CFL gap to be Aq = 25 McV, mak- 
ing A2SC = 2^/^Ao = 31.5 MeV. Recalling that 
Sfj, = M'/(8yu), this puts the second order phase 
transition at 



Ml 



6.03 A2SC = 7.60 Ao = 190.0 MeV 



(76) 

(The authors of Refs. US 0| neglected to notice 
that it is A2SC, rather than the CFL gap Aq, that 
occurs in Eqs. (|30|l and (|48|) and therefore controls 
the Sfi at which a = 0. In analyzing the crystalline 
phase in isolation, this is immaterial since either Aq 
or A2SC could be taken as the parameter defining the 
strength of the interaction between quarks. How- 
ever, in Section VI. E we shall compare the free en- 
ergies of the CFL, gCFL and crystalline phases, and 
in making this comparison it is important to take 
into account that A2SC — 2^/'^Ao.) For values of 
M"^ / fjL that are smaller than (|76() (that is, lower den- 
sities), a < and the free energy is minimized by a 
nonzero A = Amin given by 



A„ 



8Pa(J^)7off 



(77) 

and thus describes a crystalline color superconduct- 
ing phase. 

If ^off < and 7eff > 0, then the free energy 
H74(l describes a first order phase transition between 
unpaired and crystalline quark matter occurring at 



32P7eff 



(78) 



At this positive value of a, the function Vl{A) has 
a minimum at A = with $7 = 0, initially rises 
quadratically with increasing A, and is then turned 
back downward by the negative quartic term before 
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being turned back upwards again by the positive sex- 
tic term, yielding a second minimum at 



6iJ.\ 



'3|^e: 



47off 



(79) 



also with $7 = 0, which describes a crystalline color 
superconducting phase. For a < a*, the crys- 
talline phase is favored over unpaired quark matter. 
Eq. (|48|1 must be used to determine the value oi Sfi, 
and hence M'^/fi, at which a = a* and the first 
order phase transition occurs. If a* <C 1, the transi- 
tion occurs at a value of / fj, that is greater than 
(|75|l by a factor (1 -|- a*). See Fig.Ofor an exphcit 
example of plots of £7 versus A for various values of 
a for one of the crystal structures that we analyze 
in Section VI. D which turns out to have a first order 
phase transition. 

A necessary condition for the Ginzburg-Landau 
approximation to be quantitatively reliable is that 
the sextic term in the free energy is small in magni- 
tude compared to the quartic, meaning that < 
S^J^'^\|3cs/%s\■ If the transition between the unpaired 
and crystalline phases is second order, then this 
condition is satisfied close enough to the transition 
where A ^ 0. However, if ^eff < and jcS > 0, 
making the transition first order, we see from (|79|l 
that at the first order transition itself A is large 
enough to make the quantitative application of the 
Ginzburg-Landau approximation marginal. This is a 
familiar result, coming about whenever a Ginzburg- 
Landau approximation predicts a first order phase 
transition because at the first order phase transition 
the quartic and sextic terms are balanced against 
each other. Even though it is quite a different 
problem, it is worth recalling the Ginzburg-Landau 
analysis of the crystallization of a solid from a liq- 
uid [50j. There too, a Ginzburg-Landau analysis pre- 
dicts a first-order phase transition and thus predicts 
its own quantitative downfall. However, it remains 
important as a qualitative guide: it predicts a body- 
centered cubic crystal structure, and most elemen- 
tary solids are body-centered cubic near their melt- 
ing point. We shall find that our Ginzburg-Landau 
analysis predicts a first order phase transition; know- 
ing that it is therefore at the edge of its quantitative 
reliability, we shall focus in Sections VI. E and VII 
on qualitative conclusions. 

If 7eff < 0, then the Ginzburg-Landau expan- 
sion of the free energy to sextic order in (|74|l is not 
bounded from below. The transition must be first 



order, with higher-than-sextic order terms making 
the free energy bounded. In this circumstance, all 
we learn from H74(l is that the transition is first order; 
we cannot obtain an estimate of the transition point 
or of A at the first order transition. Even though 
7 is negative for many crystal structures [s^l ] , in all 
the three-fiavor crystalline phases that we present in 
Section VI. D we find that 7322 is positive and suffi- 
ciently large that 70^ = 27-1- 27322 is positive. We 
therefore need not discuss the %s < case any fur- 
ther. 



B. Two plane wave structure 

We begin with the simplest three-flavor "crystal" 
structure in which {q2} and {q3} each contain only 
a single vector, yielding a condensate 



e2^i="-A2e2a;3e2y +e2*'»^'""A3e3„/3e3„ , (80) 



in which the (us) and (ud) condensates are each 
plane waves. As explained in the previous subsec- 
tion, we shall seek solutions with A2 = A3 = A. 
We begin with such a simple ansatz both because it 
has been analyzed previously in Refs. [s^ and 
because it will yield a qualitative lesson which will 
prove extremely helpful in winnowing the space of 
multiple plane wave crystal structures. 

Let us now walk through the evaluation of all the 
coefficients in the free energy (|74f) for this two-plane 
wave structure. First, P = I (one vector in each 
of {q2} and {q3}) and as always a{Sfi) is given by 
H48|) . Next, we obtain the results for ^2 = and 
72 = 73 from the analysis of the single plane wave 
condensate in the two flavor model of Ref . [33 : 



72 



4rf -1 
1 77^ + 3 



32 (772 _ 1)3 



= 0.569 

1.637 . 



(81) 



We now turn to ^32 and 7322 which describe the in- 
teraction between the (us) and (ud) condensates and 
which we have calculated in Section V. In general, 
;532 is given by (|5^ but in this instance since {q2} 
and {qsl each contain only a single vector the sum 
in this equation reduces simply to 



P32 = ^32(0) 



(82) 



where is the angle between q2 and qg and where 
'^32(0) is given in Eq. (|58|l . /332 is plotted as a func- 
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FIG. 3: P32{4>) ~ Js2{<t>) for the two plane wave "crystal" 
structure with condensate 181)11 . is the angle between 
q2 and qa. For more complicated crystal structures, {3^2 
is given by the sum in 1591 . meaning that it is a sum of 
J32 (0) evaluated at various values of (f> corresponding to 
the various angles between a vector in {q2} and a vector 
in {qs}- 



tion of in Fig. |2| For this simple crystal structure, 
j3^2 was calculated previously in Refs. ^3iE|- 

Turning to 7322, this is given by 



7322 = --f!:322(q3,q3;q2;q2>q2 7q2) (83) 



where -^^322 is given by Eq. H69|l . As occurred in the 
evaluation of (3^2 , the sum over q- vectors in the gen- 
eral expression has reduced to evaluating -^"322 
just once, because {q2} and {q3} each contain only 
a single vector. For the special case where the last 
four arguments of K322 are the same, as in 
-^322 depends only on ^, the angle between and 
qg, and the integrals in can all be evaluated 
analytically, yielding 



322(0) 



64 (7? cos I) (772 sin2|-l 
rj^ arctan(6(0)) sin^ ^ 4 



3/2 



(r^^sin^l 



1 j Ki') b{4>) {b{(t>f - 1) 



(84) 



1 + h{(j)f 



where 



rf- sin f — 1 



77 cos ■ 



(85) 



7322 is plotted as a function of (j) in Fig. 0] 

We note that for any angle 0, both (3^2 and 27322 
are positive quantities which when added to the pos- 
itive 2/3 and 27 give positive /3off and 7eff, respec- 
tively. Hence, we see that upon making this two 
plane wave "crystal" structure ansatz we find a sec- 
ond order phase transition between the crystalline 
and unpaired phases, for all choices of the angle 0. 
We also note that both ^832(0) and 7322(0) increase 
monotonically with 0, and diverge as — *■ tt. This 
tells us that within this two plane wave ansatz, the 
most favorable orientation is = 0, namely (\2 \\ q3. 



Making this choice yields the smallest possible PcS 
and 7off within this ansatz, and hence the largest 
possible A and condensation energy, again within 
this ansatz. The divergence at — > tt tells us 
that choosing qj and q3 precisely antiparallel ex- 
acts an infinite free energy price in the combined 
Ginzburg-Landau and weak-coupling limit in which 
A ^ 5^,, Aq <C /i, meaning that in this limit if 
we chose = tt we find A = 0. Away from the 
Ginzburg-Landau limit, when the pairing rings on 
the Fermi surfaces widen into bands, choosing = tt 
exacts a finite price meaning that A is nonzero but 
smaller than that for any other choice of 0. All these 
results confirm conclusions drawn in Refs. [s^ ^| 
based only upon the result for /332(0). 

The high cost of choosing qj and qg precisely an- 
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FIG. 4; 7322(0) for the two plane wave "crystal" struc- 
ture with condensate (I8UII . 4> is the angle between q2 and 
qa- 7322(0) = 0.243 and 7322 (<?!') increases monotonically 
with (f>. 



tiparallel can be understood qualitatively as arising 
from the fact that in this case the ring of states 
on the w-quark Fermi surface that "want to" pair 
with d-quarks coincides precisely with the ring that 
"wants to" pair with s-quarks [41|. (For example, 
if (\2 and q3 point in the —z and +z directions, A2 
(A3) describes pairing between s-quarks (d-quarks) 
within a ring in the northern hemisphere of the s- 
((i-)Fermi surface and w-quarks within a ring in the 
southern hemisphere of the u-Fermi surface. The 
rings on the M-Fermi surface coincide, as illustrated 
in Fig. 2 of Ref. ^J.) In the most favorable case 
within the two-plane wave ansatz, where || 
the two pairing rings on the M-quark Fermi surface 
are centered on antipodal points (^^^^ exam- 

ple, if c\2 and q3 both point in the direction, A2 
(A3) describes pairing of s-quarks (d-quarks) within 
a ring in the southern (northern) hemisphere of the 
s- (d-)Fcrmi surface and u-quarks within rings in the 
(northern) southern hemisphere of the w-Fermi sur- 
face.) 

The simple two plane wave ansatz (I8UII has been 
analyzed in the same NJL model that we employ 
upon making the weak-coupling approximation but 
without making a Ginzburg-Landau approximation 
in Ref. the qualitative lessons that we 

have learned from the Ginzburg-Landau approxima- 
tion remain valid, including the favorability of the 
choice = 0, but we learn further that in the two 
plane wave case the Ginzburg-Landau approxima- 



tion always underestimates A |4l| . We also see from 
Ref. li^l that the A at which the Ginzburg-Landau 
approximation breaks down shrinks as </> — > tt. We 
can understand this result as follows. The sextic 
term in the free energy l|7Hl is small compared to 
the quartic term only if A^ ^ 5/x^/3eff/7oft, making 
this a necessary condition for the quantitative va- 
lidity of the Ginzburg-Landau approximation. As 
4> —I- TT, 7ofr diverges more strongly than /?off: from 
and we find that as ^ tt, 

3- 3^ ( ^ y 

7cff 322 1^^^^^^^ J . 

(86) 

Therefore the Ginzburg-Landau calculation predicts 
that its own breakdown will occur at a A that de- 
creases with increasing (f>, as found in Ref. by 
explicit comparsion with a calculation that does not 
employ the Ginzburg-Landau approximation. 

C. Implications for more plane waves: 
qualitative principles for favorable crystal 
structures 

In this subsection we ask what lessons we can learn 
from the evaluation of the Ginzburg-Landau coeffi- 
cients for the two plane wave "crystal" structure in 
Section VLB for crystal structures with more than 
one vector in {q2} and {q3}. 

First, we can conclude that P32 is positive for any 
choice of {q2} and {q3}. The argument is simple: 
P32 is given in general by H59|l . a sum over J32 eval- 
uated at a host of angles corresponding to all angles 
between a vector in {q2} and {q^^}. But, we see from 
Fig-Elthat J32 is positive at any angle. 

Second, we cannot draw such a conclusion about 
7322. This coefficient is a sum over contributions 

of the form i?322(qt Q^, q2, Q^, Q^) where the 
last four momentum vector arguments, selected from 
{q2}, must satisfy (|62|l . The calculation in Section 
VLB whose result is plotted in Fig. 0] only demon- 
strates that those contributions in which the four 
q3 arguments are selected to all be the same vec- 
tor are positive. For any crystal structure in which 
{q2} contains two or more vectors, there are other 
contributions to 7322 that we have not evaluated in 
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this section which depend on one vector and sev- 
eral Qg vectors, and thus on more than one angle. 
We know of instances where individual contributions 
-^322(q3, Qa, *l2 ' Q2) Crystal structures that 

we describe below are negative. However, we have 
found no crystal structure for which 7322 is negative. 

The final lesson wc learn is that crystal structures 
in which any of the vectors in {q2} arc close to an- 
tiparallel to any of the vectors in {q^} are strongly 
disfavored. (The closer to antiparallel, the worse, 
with the free energy penalty for A 7^ diverging for 
the precisely antiparallel case, driving A to zero.) If 
a vector in {q2} is antiparallel (or close to antipar- 
allel) to one in {qg}, this yields infinite (or merely 
large) positive contributions to P32 and to 7322 and 
hence to peS and ^^ff- In the case of P32, these large 
positive contributions cannot be cancelled since all 
contributions are positive. In the case of 7322, neg- 
ative contributions are possible but we know of no 
instances of divergent negative contributions to 7322 
or indeed to any other coefficient in the Ginzburg- 
Landau expansion. The divergent positive contribu- 
tions are associated with the tangential intersection 
(in the case of /3 and 7 js^) or coincidence (in the 
case of of /?32 and 7322) of pairing rings on Fermi 
surfaces. We know of no configuration of rings that 
leads to an infinitely favorable (as opposed to un- 
favorable) free energy in the combined Ginzburg- 
Landau and weak-coupling limits. So, although we 
do not have a proof that the divergent positive con- 
tributions to 7322 arising as vectors in {q2} and {q3} 
approach one another's antipodes are uncancelled, 
we also see no physical argument for how this could 
conceivably arise. Certainly in all example crystal 
structures that we have considered, /332 and 7322 
and hence f3eS and 7off diverge as vectors in {q2} 
and {q3} approach one another's antipodes. 

We can now summarize the qualitative principles 
that we have arrived at for constructing favorable 
crystal structures for three-flavor crystalline color 
superconductivity. First, as described in Section 
III the sets {q2} and {q3} should each separately 
be chosen to yield crystal structures which, seen as 
separate two-flavor crystalline phases, are as favor- 
able as possible. In Section III we have reviewed 
the results of Ref. [s^ for how this should be done, 
and the conclusion that the most favored {q2} or 
{q3} in isolation consists of eight vectors pointing 
at the corners of a cube. Second, the new addition 
in the three-flavor case is the qualitative principle 
that {q2} and {q3} should be rotated with respect 



to each other in such a way as to best keep vectors 
in one set away from the antipodes of vectors in the 
other set. 



D. Multiple plane waves 

In Table I we describe 11 different crystal struc- 
tures that we have analyzed, and in Table II we give 
the coefficients that specify each Ginzburg-Landau 
free energy ((7^ . T he B 's and 7's were calculated 
as described in Ref. [33; /332's and 7322 's were cal- 
culated as described in Section V. We also give the 
combinations /3cff and 7off defined in H74|l that spec- 
ify the free energy as in H75|l . In those cases in which 
PcB < 0, the phase transition between the crystalline 
phase and the unpaired phase is first order, occur- 
ring where a = a* with a* given by (|78|l . At the 
first order phase transition, unpaired quark matter 
with A = and crystalline quark matter with A (a*) 
given in (|79|l have the same free energy. We give 
both a, and A(a*) in Table II. 

The first row of the Tables describes the sim- 
ple "crystal structure" analyzed in detail in Section 
VLB, in which both {q2} and {q3} contain just a 
single vector, with q2 || q3 as we have seen that 
this is the most favorable choice for the angle be- 
tween q2 and qg. This condensate carries a baryon 
number current which means that the unpaired gap- 
less fermions (in "blocking regions" in momentum 
space 113, 1^ ) must carry a current that is equal in 
magnitude but opposite in direction |30l | . The anal- 
ysis of this "crystal structure" in Sections VLB and 
VI.C has proved instructive, giving us qualitative in- 
sight that we shall use to understand all the other 
crystal structures. However, in all rows in the Tables 
other than the first we have chosen crystal structures 
with condensates that carry no net current, mean- 
ing that the gapless fermions need carry no current. 
There is nothing in our mean-field analysis that pre- 
cludes condensates carrying a net current, but we 
do not analyze them here primarily because it sim- 
plifies our task but also because we expect that, be- 
yond mean-field theory, a phase containing gapless 
fermions carrying a net current is unlikely to be the 
favored ground state. 

Let us next examine the last two rows of the 
Tables. Here, we consider two crystal structures 
in which {q2} and {q3} each contain eight vectors 
forming cubes. Since the cube is the most favorable 
two-fiavor crystal structure according to the analy- 
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Structure 


Description 


Largest Angle 


2PW 


{q2} and {qs} coincide; each contains one vector. (So, 2 plane waves with q2 | qa.) 


0° 


SqX 


{q2} and {qa} each contain two antiparallel vectors. The four vectors togetlier form 
a square; those from {q2} and those from {qa} each form one stroke of an "X". 


90° 


Tetrahedron 


{q2} and {qg} eacli contain two vectors. The four together form a tetrahedron. 


109.5° 


2Triangles 


{q2} and {q^} coincide; each contains three vectors forming a triangle. 


120° 


Cube X 
See Eq. (TO 


jool and Iq-il each contain 4 vectors forming a rectangle. The 8 vectors together 
form a cube. The 2 rectangles intersect to look like an "X" if viewed end-on. 


109.5° 


2Tet 


{q2} and {q^} coincide; each contains four vectors forming a tetrahedron. 


109.5° 


Twisted 
Cube 


{q2} and {qs} each contain four vectors forming a square which could be one face of 
a cube. Instead, the eight vectors together form the polyhedron obtained by twisting 
the top face of a cube by 45° relative to its bottom face. 


143.6° 


2Octa90xy 


jqo} and Iqij each contain 6 vectors forming an octahedron. The |qo| vectors point 
along the positive and negative axes. The {qsj-octahedron is rotated relative to the 
{q2}-octahedron by 90° about the (1, l,0)-axis. 


135° 


20cta45xyz 


{q2} and {qa} each contain 6 vectors forming an octahedron. The {q2} vectors point 
along the positive and negative axes. The {qaj-octahedron is rotated relative to the 
{q2}-octahedron by 45° about the (1, 1, l)-axis. 


143.6° 


2Cube45z 
See Eq. (|HZJ 


{q2} and {qs} each contain 8 vectors forming a cube. The {q2} vectors point along 
(±1,±1,±1). The {qal-cube is rotated relative to that by 45° about the z-axis. 


143.6° 


2Cube45xy 


{q2} and {qs} each contain 8 vectors forming a cube. The {q2} vectors point along 
(±1,±1,±1). The {qal-cube is rotated relative to that by 45° about the (1, l,0)-axis. 


154.5° 



TABLE L Descriptions of the crystal structures whose Ginzburg-Landau coefficients are given in Table II. The third 
column is the largest angle between any vector in {q2} and any vector in {qs}. Other things being equal, we expect 
that the larger the largest angle, meaning the closer vector(s) in {q2} get to vector(s) in {qs}, the bigger the P32 and 
7322 and hence the bigger the j3cs and 7cff , and hence the less favorable the structure. 



Structure 




^32 




7 


7322 


7cff 


Of* 


A(a,)/5M 


2PW 


0.569 


0.250 


1.388 


1.637 


0.243 


3.760 








SqX 


0.138 


1.629 


1.906 


1.952 


2.66 


9.22 








Tetrahedron 


-0.196 


2.146 


1.755 


1.450 


7.21 


17.29 








2Triangles 


-1.976 


4.647 


0.696 


1.687 


13.21 


29.80 








CubcX 


-10.981 


6.961 


-15.001 


-1.018 


19.90 


37.76 


0.140 


0.548 


2Tet 


-5.727 


7.439 


-4.015 


4.350 


30.35 


69.40 


0.0054 


0.208 


Twisted Cube 


-16.271 


12.445 


-20.096 


-37.085 


315.5 


556.8 


0.0170 


0.165 


2Octa90xy 


-31.466 


18.665 


-44.269 


19.711 


276.9 


593.2 


0.0516 


0.237 


20cta45xyz 


-31.466 


19.651 


-43.282 


19.711 


297.7 


634.9 


0.0461 


0.226 


2Cube45z 


-110.757 


36.413 


-185.101 


-459.24 


1106. 


1294. 


0.310 


0.328 


2Cube45xy 


-110.757 


35.904 


-185.609 


-459.24 


11358. 


21798. 


0.0185 


0.0799 



TABLE II: Ginzburg-Landau coefficients for three-flavor crystalline color superconducting phases with various crystal 
structures, described in Table I. q, is the a at which the transition from unpaired quark matter to a given crystalline 
phase occurs: q, = if (3cH > and the transition is second order; a, is given by (1781 if /3off < and the transition 
is first order. For a first order transition, A(q,), given in l|79|l . is the magnitude of the gap at the transition. 
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sis of Ref. , evident in the large negative P and 
7 for both these crystal structures in Table II, this 
should be a good starting point. We cannot have the 
two cubes coincident, as in that case there are vec- 
tors from {q2} and vectors from {qs} separated by 
a 180° angle, yielding infinite positive contributions 
to both P32 and 7322- So, we rotate the {qgj-cube 
relative to the {q2} cube, in two different ways in 
the 2Cubc45z and 2Cubc45xy crystal structures de- 
scribed in Table I. 

We explain explicitly in Appendix B why translat- 
ing one cube relative to the other in position space 
by a fraction of a lattice spacing does not allevi- 
ate the problem: a relative rotation of the {us) and 
(ud) condensates is required. Qualitatively, this re- 
flects the nature of the difficulty that occurs when 
a {q2} vector is opposite to a {Qs} vector. It can 
be thought of as arising because the (us) and (ud) 
condensates both want to "use" those up quarks ly- 
ing on the same ring on the up Fermi surface. It 
therefore makes sense that a relative rotation is re- 
quired. Quantitatively, what we show in Appendix 
B is that n does not change if we translate the {us) 
condensate relative to the (ud) condensate. 

We see in Table I that in the 2Cube45z structure, 
the largest angle between vectors in {q2} and {q3} 
is 143.6° whereas in the 2Cube45xy structure, that 
largest angle is 154.5° meaning that the rotation we 
have chosen does a less good job of keeping {q2}- 
vectors away from the antipodes of {q3} vectors. 
Correspondingly, we see in Table II that 2Cube45xy 
has a much larger 7322 and hence 7off , and hence has 
a first order phase transition occurring at a smaller 
a* and with a smaller A (a*). This is an example 
confirming our general principle that, other things 
being equal, crystal structures in which {q3} vectors 
come closer to {q2} vectors will be disfavored. Ac- 
cording to this principle, the 2Cube45z crystal struc- 
ture should be particularly favorable as it employs 
the relative rotation between the two cubes that does 
the best possible job of keeping them apart. 

We now turn to crystal structures with fewer than 
16 plane waves. By having fewer than 8 plane waves 
in {q2} and {q3}, we arc no longer optimizing the 
two-flavor /3 and 7. However, with fewer vectors it is 
possible to keep the {q2}- and {qgj-vectors farther 
away from each other's antipodes. We list two crys- 
tal structures in which {q2} and {qg} have 6 waves 
forming octahedra. These are not particularly favor- 
able two-flavor structures — 7 is positive rather than 
being large and negative for the cube. 20cta45xyz 



has the same largest angle between {q2}- and {qs}- 
vectors as 2Cube45z, but its significantly more pos- 
itive P and 7 make it significantly less favorable. 
Choosing the 2Octa90xy structure instead reduces 
the largest angle between {q2}- and {q3}-vectors 
from 143.6° to 135°, which improves PcS and 7off, 
but only slightly. 

We investigate three crystal structures in which 
{q2} and {q3} each contain 4 vectors. Among these, 
the Twisted Cube is strongly disfavored by its signif- 
icantly larger largest angle between {q2}- and {qg}- 
vectors. CubeX and 2Tet are both constructed by 
choosing {q2} and {qs} as subsets containing half 
the vectors from a cube. In the 2Tet structure, 
we choose the tetrahedra coincident since this does 
the best job of keeping vectors in {q2} and {q3} 
away from each other's antipodes. (Choosing the 
two tetrahedra so that their union forms a cube is 
the worst possible choice, as vectors in {q2} and 
{q3} are then antipodal.) In the CubeX structure, 
we choose the two rectangles such that their union 
forms a cube, as this does the best job of reduc- 
ing the largest angle between vectors in {q2} and 
{q3}; making the rectangles coincident would have 
been the worst possible choice. CubeX and 2Tet 
have the same largest angle, but they differ consider- 
ably in that the {q2} and {q^} rectangles that make 
up CubeX are more favorable two-flavor structures 
(lower P and 7) than the tetrahedra that make up 
2Tet. We see from Table II that the CubeX struc- 
ture, with only 8 vectors in total, is particularly 
favorable: it is not possible to tell from Table II 
whether it is more or less favorable than 2Cube45z, 
since one has the larger a* while the other has the 
larger A(a*). We shall evaluate their free energies 
below, and confirm that they are indeed compara- 
ble, and that these two structures have the lowest 
free energy of any in the Tables. 

In the remaining rows of the Tables, we investigate 
one crystal structure in which {q2} and {q3} each 
contain 3 vectors, and two in which each contain 
2 vectors. These structures all have positive /3off, 
and hence second order phase transitions, and so 
are certainly not favored. 

Inspecting the results in Table II shows that in 
all cases where we have investigated different three- 
flavor crystal structures built from the same {q2} 
and {q3}, the one with the relative rotation be- 
tween the two polyhedra that yields the smaller 
largest angle between vectors in {q2} and {q3} is 
favored. And, in all cases where we have investi- 
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gated two crystal structures with same largest an- 
gle between vectors in {q2} and {q3}, the one built 
from the more favorable two-flavor crystal structure 
is favored. We thus find no exceptions to the quali- 
tative principles we described in Section VI. C. How- 
ever, these qualitative principles certainly do not 
explain all the features of the results in Table II. 
For example, we have no qualitative understanding 
of why 2Cube45z and 2Cubc45xy have such simi- 
lar /332, whereas 2Cube45xy has a much larger 7322 
as expected. For example, we have no qualitative 
understanding of why 7322 increases much more in 
going from 2Cube45xy to 2Cube45z than it does in 
going from 2Octa90xy to 20cta45xyz. The calcula- 
tions must be done; the qualitative principles are a 
good guide, but not a substitute. 

The final crystal structure that we describe is one 
in which {q2} is a cube while {q3} is an octahedron, 
with the six {q3}-vectors pointing at the centers of 
the faces of the {q2}-cube. So, if the {qjj-vectors 
are taken to point along the (±1, ±1, ±1) directions 
then the {q3} vectors point along the positive and 
negative axes. We chose to investigate this struc- 
ture because it seems particularly symmetric and 
because it has an unusually small largest angle be- 
tween vectors in {q2} and {q3} given the large num- 
ber of vectors in total: 125.3°. Because {q2} and 
{q3} are not congruent, 7^ P3 and 72 7^ 73. AH 
these coefficients can be found in Table II. Wc find 
^32 = 24.5 1 , 7322 = 4 1 9.9 and 7233 = 4943. Be- 
cause {q2} and {93} are not exchange symmetric, 
the general argument that we gave in Section VI. A 
for why extrema of ri(A2, A3) — i.e. solutions to the 
gap equations — occur at A2 = A3 does not apply. 
However, we find that at the solution A2 and A3 
differ by less than 20%. The large values of 7233 and 
7322 make this crystal structure quite unfavorable — 
even though it has a (weak) first order phase transi- 
tion, its free energy turns out to be comparable only 
to that of the 2PW structure, far above the free en- 
ergy of the favored CubeX and 2Cube45z structures. 
Furthermore, the arguments of Appendix A do not 
apply to a crystal structure like this, meaning that 
we do not expect this solution with A2 7^ A3 to be 
neutral. For this reason, and because it appears to 
be free-energetically unfavorable anyway, we will not 
investigate it further. We cannot say whether choos- 
ing {q2} and {q3} to not be exchange symmetric 
generically yields an unfavorable crystal structure, 
as we have not investigated many possibilities. 

We have certainly not done an exhaustive search 



of three-flavor crystal structures. For example, 
we have only scratched the surface in investigat- 
ing structures in which {q2} and {q3} are not ex- 
change symmetric. We have investigated the struc- 
tures that are the best that we can think of accord- 
ing to the qualitative principles described in Section 
VI. C. Readers should feel free to try others. (We are 
confident that in 2Cubc45z wc have found the most 
favorable structure obtained by rotating one cube 
relative to another. We arc not as confident that 
CubeX is the best possible structure with fewer than 
8-1-8 vectors.) As we shall see in Section VI. E, how- 
ever, the two most favorable structures that we have 
found, 2Cube45z and CubeX, are impressively ro- 
bust and do a very good job of making the case that 
three-flavor crystalline color superconducting phases 
are the ground state of cold quark matter over a wide 
range of densities. If even better crystal structures 
can be found, this will only further strengthen this 
case. 



E. Free energy comparisons 

We can now evaluate and plot the gap parameter 
A and free energy n{A) for all the crystal structures 
described in Table I, whose Ginzburg-Landau coeffi- 
cients are given in Table II. For a given crystal struc- 
ture, 0(A) is given by Eq. ifT^ . with (3cS and 7eff 
taken from Table II. The quadratic coefficient a is 
related to Sfi by Eq. (gHJ). Recall that we have made 
the approximation that 6^2 = = S/i = / (8/i), 
valid up to corrections of order / . At any value 
of Mj/ji, we can evaluate a((5/i) and hence 0(A), 
determine A by minimizing O, and finally evaluate 
the free energy O at the minimum. In Fig. [SJ we 
give an example of 0(A) for various M^/jj, for one 
crystal structure with a first order phase transition 
(CubcX), illustrating how the first order phase tran- 
sition is found, and how the A solving the gap equa- 
tions — i.e. minimizing O — is found . We plot A 
and at the minimum versus M"^ / ji in Figs. El and 
[3 for some of the crystal structures in Tables I and 
II. 

In Figs. El and [7| we show two examples of crys- 
tal structures for which the phase transition to the 
unpaired state is second order: 2PW and SqX. (See 
Table I for descriptions of these structures.) The 
second order phase transition occurs at M'^/fi = 
7.6OA0 = 190.0 McV, where a = 0. (See Eq.' lfTC|) .') 
We show four examples of crystal structures with 
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FIG. 5: Free energy f2 vs. A for the CubeX crystal structure, described in Table II, at four values of /ji. From top 
curve to bottom curve, as judged from the left half of the figure, the curves are //i = 240, 218.61, 190, and 120 MeV, 
corresponding to a = 0.233, 0.140, 0, -0.460. The first order phase transition occurs at M^/^i = 218.61 MeV. The 
values of A and Q at the minima of curves like these are what we plot in Figs. |S| and [7| 



first order phase transitions, occurring where a — 
a* > meaning at some / fj, > 190.0 MeV. We 
show the two most favorable structures that we have 
found: CubeX and 2Cubc45z. And, we show two ex- 
amples (2Tet and 2Octa90xy) of structures with first 
order phase transitions that are more favorable than 
the structures with a second order transition, but 
less favorable than CubeX and 2Cubc45z. 

In Figs. El and we have chosen the interaction 
strength between quarks such that the CFL gap pa- 
rameter at Afg = is Ao = 25 MeV. However, our 
results for both the gap parameters and the con- 
densation energy for any of the crystalline phases 
can easily be scaled to any value of Aq. We saw 
in Section V that the quartic and sextic coefficients 
in the Ginzburg-Landau free energy do not depend 
on Ao. And, recall from Eq. H48(l that Aq enters 
a only through the combination A2sc/<^Aij where 
A2SC = 23 Ao and dfi = M^/{8fi). This means that 
if we pick a Ao ^ 25 MeV, the curves describing the 
gap parameters for the crystalline phases in Fig. 
are precisely unchanged if we rescale both the verti- 



cal and horizontal axes proportional to Ao/25 MeV. 
In the case of Fig. |3 the vertical axis must be 
rescaled by (Ao/25 McV)^. Of course, the weak- 
coupling approximation Ao fj., which we have 
used for example in simplifying the propagators in 
(|52() . will break down if we scale Ao to be too large. 
Wc cannot evaluate up to what Ao we can scale 
our results reliably without doing a calculation that 
goes beyond the weak-coupling limit. However, such 
calculations have been done for the gCFL phase in 
Ref. 23 1 where it turns out that the gaps and con- 
densation energies plotted Figs. El and [3 scale with 
Ao and Aq to good accuracy for Ao < 40 MeV with 
/i = 500 MeV, but the scaling is significantly less ac- 
curate for Aq = 100 McV. Of course, for Aq as large 
as 100 McV, any quark matter in a compact star is 
likely to be in the CFL phase. Less symmetrically 
paired quark matter, which our results suggest is in 
a crystalline color superconducting phase, will occur 
in compact stars only if Ao is smaller, in the range 
where our results can be expected to scale well. 

The qualitative behavior of A at smaller M^//x, 
well to the left of the unpaired/crystalline phase 
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FIG. 6: Gap parameter A versus /fi for three-flavor crystalline color superconducting phases with various crystal 
structures. The crystal structures are described in Table II. For comparison, we also show the CFL gap parameter and 
the gCFL gap parameters Ai, A2 and A3 [TH. Il2|. Recall that the splitting between Fermi surfaces is proportional 
to M^/fi, and that small (large) Atg/fJ, corresponds to high (low) density. 



transitions in Fig. |S1 can easily be understood. The 
quadratic, quartic and sextic coefficients in the free 
energy (jTU are a(i5/x), f3es = Pcs/SfJ-'^ and 7eff = 
7eff/^M^- If tended to a constant at small 6fi, then 
the solution A,„i,i that minimizes would be pro- 
portional to (5/Lt. (Sec Eq. Gil).) In fact, from 
we see that a oc log Sfi at small 5/^, meaning that, 
according to lf77|l . Amin should vanish slightly more 
slowly than linear as Mj//i oc (5/i ^ 0, as in Fig.|H| 
And, since the A's vanish for Sfi 0, so do the 
condensation energies of Fig. [7| 

Fig. |B| can be used to evaluate the validity of 
the Ginzburg-Landau approximation. The simplest 
criterion is to compare the A's for the crystalline 
phases to the CFL gap parameter Ao. This is the 
correct criterion in the vicinity of the 2nd order 
phase transition point, where Sfi = M^/(8/i) ~ Aq. 
Well to the left, it is more appropriate to compare 
the A's for the crystaUine phase to Sfi = M^/{8^). 
By either criterion, we sec that all the crystal struc- 
tures with first order phase transitions (including the 



two that are most favored) have A's that are large 
enough that the Ginzburg-Landau approximation is 
at the edge of its domain of validity, a result which 
wc expected based on the general arguments in Sec- 
tion VI.A. Note that the Ginzburg-Landau approxi- 
mation is controlled for those structures with second 
order phase transitions only near the second order 
phase transition, again a result that can be argued 
for on general grounds. 

Fig. [71 makes manifest one of the central conclu- 
sions of our work. The three-flavor crystalline color 
superconducting phases with the two most favored 
crystal structures that we have found are robust 
by any measure. Their condensation energies reach 
about half that of the CFL phase at Ms = 0, remark- 
able given that in the CFL phase pairing occurs over 
the whole of all three Fermi surfaces. Correspond- 
ingly, these two crystal structures are favored over 
the wide range of seen in Fig. [T] and given in 

Eq. @. 

Taken literally. Fig. indicates that within the 
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FIG. 7: Free energy Q, versus /fi for the three-flavor crystaUine color superconducting phases with various crystal 
structures whose gap parameters are plotted in Fig. 15] The crystal structures are described in Table II. Recall that 
the gCFL phase is known to be unstable, meaning that in the regime where the gCFL phase free energy is plotted, 
the true ground state of three-flavor quark matter must be some phase whose free energy lies below the dashed line. 
We see that the three-flavor crystalline color superconducting quark matter phases with the most favorable crystal 
structures that we have found, namely 2Cube45z and CubeX described in 18711 and 19011 . have sufficiently robust 
condensation energy (sufficiently negative Q) that they are candidates to be the ground state of three-flavor quark 
matter over a wide swath of Mf//i, meaning over a wide range of densities. 



regime © of the phase diagram occupied by crys- 
talline color superconducting quark matter, the 
2Cube45z phase is favored at lower densities and 
the CubeX phase is favored at higher densities. Al- 
though, as detailed in Sections VI.C and VI. D, we 
do have qualitative arguments why 2Cube45z and 
CubeX are favored over other phases, we have no 
qualitative argument why one should be favored over 
the other. And, we do not trust that the Ginzburg- 
Landau approximation is sufficiently quantitatively 
reliable to trust the conclusion that one phase is fa- 
vored at higher densities while the other is favored 
at lower ones. We would rather leave the reader with 
the conclusion that these are the two most favorable 
phases we have found, that both are robust, that 
the crystalline color superconducting phase of three- 
flavor quark matter with one crystal structure or the 
other occupies a wide swath of the QCD phase dia- 



gram, and that their free energies are similar enough 
to each other that it will take a beyond-Ginzburg- 
Landau calculation to compare them reliably. 



VII. CONCLUSIONS, IMPLICATIONS, AND 
FUTURE WORK 



We have evaluated the gap parameter and free 
energy for three-flavor quark matter in crystalline 
color superconducting phases with varied crystal 
structures, within a Ginzburg-Landau approxima- 
tion. Our central results are shown in Figs.|Hland[3 
Descriptions of the crystal structures that we have 
investigated, together with the coefficients for the 
Ginzburg-Landau free energy I74|) for each struc- 
ture, are given in Tables I and II. 

We have found two qualitative rules that guide our 
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understanding of what crystal structures are favored 
in three-flavor crystalhne quark matter. First, the 
{ud) and {us) condensates separately should be cho- 
sen to have favorable free energies, as evaluated in 
the two-flavor model of Ref. [33 • Second, the {ud) 
and {us) condensates should be rotated relative to 
each other in such a way as to maximize the an- 
gles between the wave vectors describing the crystal 
structure of the {ud) condensate and the antipodes 
of the wave vectors describing the {us) condensate. 
This second qualitative rule can be understood as 
minimizing the "competition" between the two con- 
densates for up quarks on the up Fermi surface, as 
first elucidated in a simpler setting in Ref. • 

Fig.dshows that over most of the range of M'^/fi 
where it was once considered a possibility, the gCFL 
phase can be replaced by a much more favorable 
three-flavor crystalline color superconducting phase. 
However, Fig. [3 also indicates that it is hard to find 
a crystal structure which yields a crystalline phase 
that has lower free energy than the gCFL phase at 
the lowest values of /fi (highest densities) in the 
"gCFL window" , closest to the CFL^gCFL transi- 
tion. This narrow window where the gCFL curve re- 
mains the lowest curve in Fig.Qis therefore the most 
likely place in the QCD phase diagram in which to 
find the gCFL phase augmented by current-carrying 
meson condensates described in Refs. ^^^3- Ex- 
cept within this window, the crystalline color su- 
perconducting phases with either the CubeX or the 
2Cube45z crystal structure provide an attractive res- 
olution to the instability of the gCFL phase. 

The three-flavor crystalline color superconducting 
phases with the CubeX and 2Cube45z crystal struc- 
tures have condensation energies that can be as large 
as half that of the CFL phase. This robustness 
makes them the lowest free energy phase that we 
know of, and hence a candidate for the ground state 
of QCD, over a wide range of densities. To give a 



sense of the implications of the range of M'^/fi over 
which crystalline color superconductivity is favored, 
given by Eq. (0) and shown in Fig. if we sup- 
pose that Ao = 25 MeV and Ms = 250 MeV, the 
window (El) translates to 240MeV < /i < 847MeV. 
With these choices of parameters, then, the lower 
part of this range of fj, (higher part of the range of 
M'^ / fi) is certainly superseded by nuclear matter. 
And, the high end of this range extends far beyond 
the fi ^ 500 MeV characteristic of the quark mat- 
ter at the densities expected at the very center of 
compact stars. Our result therefore suggests that if 
compact stars have quark matter cores, it is entirely 
reasonable to suppose that the entire quark matter 
core could be in a crystalline color superconducting 
phase. Of course, if Aq is larger, say ^ 100 MeV, the 
entire quark matter core could be in the CFL phase. 
And, there are reasonable values of Aq and Ms for 
which the outer layer of a possible quark matter core 
would be in a crystalline phase while the inner core 
would not. We do not know Aq and Mg well enough 
to answer the question of what phases of quark mat- 
ter occur in compact stars. However, our results add 
the possibility that as much as all of the quark mat- 
ter in a compact star could be in a crystalline color 
superconducting phase to the menu of options that 
must ultimately be winnowed by confrontation with 
astrophysical observations. 

We have identified two particularly favorable crys- 
tal structures, using the qualitative rules described 
above and by direct calculation. We do not believe 
that our Ginzburg-Landau approximation is suffi- 
ciently accurate to trust its determination of which 
of these two structures is more favorable. For this 
reason, we wish to leave the reader with a picture 
of both the 2Cube45z and CubeX crystal structures 
in position space. In the 2Cube45z phase, the color- 
flavor and position space dependence of the conden- 
sate, defined in (|20|l and (|21|) . is given by 
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cos — ix + y + z) + cos — (— x + y + z) 
a a 
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+ cos — [x ~ y + z) + cos — [—X — y + z) 
a a 



2A 



cos — (^V2x + z^ + cos — (V2y + 



(87) 



2tt 



-\f2y + + cos — [^\f2x + z 



where a and /? (i and j) are color (flavor) indices 
and where 



a = 



VStt _ 4.536 
6fi 



1.764M? 



is the lattice spacing of the face-centered cubic 
crystal structure. For example, with M'^/fi — 
100,150,200 MeV the lattice spacing is a = 
72,48,36 fm. Eq. (|57|l can equivalently be written 
as 

AcF{x)ai,l33 = e2a/3e2yA2(r) + e^a 13^31] A3 (r) , 

(89) 

with l|57|l providing the expressions for A2(r) and 
A3(r). A three-dimensional contour plot that can be 



seen as depicting either A2(r) or A3(r) separately 
can be found in Ref. [s^. We have not found an 
informative way of depicting the entire condensate 
in a single contour plot. Note also that in H87|l and 
below in our description of the CubcX phase, we 
make an arbitrary choice for the relative position of 
A3(r) and A2(r). We show in Appendix B that one 
can be translated relative to the other at no cost 
in free energy. Of course, as we have investigated 
in detail in Section VI, rotating one relative to the 
other changes the Ginzburg-Landau coefficients ^32 
and 7322 and hence the free energy. 

In the CubeX phase, the color-flavor and position 
space dependence of the condensate is given by 



J 



AcF{x)ai,Pj =e2a/3e2ij 2A 



2tt ^ N 27r ^ 

cos — [x + y + z) + cos — (—X — y + z) 



2A 



cos — (~x + y + z) + cos — (x ^ y + z) 



(90) 



We provide a depiction of this condensate in Fig. |S1 

The gap parameter A is large enough in both the 
2Cube45z and CubeX phases that the Ginzburg- 
Landau approximation that we have used to obtain 
our results is being pushed to the limits of its valid- 
ity. Therefore, although we expect that the qualita- 
tive lessons that we have learned about the favorabil- 



ity of crystalline phases in three-flavor quark matter 
are valid, and expect that the relative favorability of 
the 2Cubc45z and CubeX structures and the qual- 
itative size of their A and condensation energy are 
trustworthy, we do not expect quantitative reliabil- 
ity of our results. There is therefore strong moti- 
vation to analyze crystalline color superconducting 
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FIG. 8: The CubeX crystal structure of Eq. The fi gure extends from to a/2 in the x, y and z directions. Both 

A2(r) and A3(r) vanish at the horizontal plane. A2(r) vanishes on the darker vertical planes, and A3(r) vanishes 
on the lighter vertical planes. On the upper (lower) dark cylinders and the lower (upper) two small corners of dark 
cylinders, A2(r) = +3. 3 A (A2(r) — — 3.3A). On the upper (lower) lighter cylinders and the lower (upper) two small 
corners of lighter cylinders, A3(r) — — 3.3A (A3(r) — +3.3A). Note that the largest value of |A/(r)| is 4A, occurring 
along lines at the centers of the cylinders. The lattice spacing is a when one takes into account the signs of the 
condensates; if one looks only at |A7(r)|, the lattice spacing is a/2, a is given in (1881 . In 191)1 and hence in this figure, 
we have made a particular choice for the relative position of A3(r) versus A2(r). We show in Appendix B that one 
can be translated relative to the other with no cost in free energy. 



quark matter with these two crystal structures with- 
out making a Ginzburg-Landau approximation. It 
will be very interesting to see whether the Ginzburg- 
Landau approximation underestimates A and the 
condensation energy for the crystalline phases with 
CubeX and 2Cube45z crystal structures, as it does 
for the much simpler 2PW structure (in which A2(r) 
and A3(r) are each single plane waves) |4l] |. 

Even prior to having a beyond-Ginzburg-Landau 
analysis available, having an ansatz (actually, two 
ansatze) for the crystal structure and a good qual- 
itative guide to the scale of A and ri(A) should 
allow significant progress toward the calculation of 
astrophysically relevant obscrvables. For example, 
it would be interesting to evaluate the effects of a 



crystalline color superconducting core on the rate 
at which a neutron star cools by neutrino emission. 
The specific heat of crystalline color superconduct- 
ing quark matter is linear with T because of the 
presence of gapless quark excitations at the bound- 
aries of the regions in momentum space where there 
arc unpaired quarks Calculating the heat ca- 

pacity of the CubeX and 2Cube45z structures should 
therefore yield only quantitative changes relative to 
that for unpaired quark matter, unlike in the gCFL 
case where the heat capacity is parametrically en- 
hanced The neutrino emissivity should turn 
out to be significantly suppressed relative to that 
in unpaired quark matter. The evaluation of the 
phase space for direct URCA neutrino emission from 
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the CubcX and 2Cubc45z phases will be a nontriv- 
ial calculation, given that thermally excited gapless 
quarks occur only on patches of the Fermi surfaces, 
separated by the (many) pairing rings. (The direct 
URCA processes u + e s + v and s ^ u + e -\- D 
require s, u and e to all be within T of a place in 
momentum space where they are gapless and at the 
same time to have Pu + Pe = Ps to within T. Here, 
T ~ keV is very small compared to all the scales 
relevant to the description of the crystalline phase 
itself.) 

Beginning with Ref. [^OIj one of the motivations 
for the study of crystalline color superconducting 
quark matter has been the possibility that, if present 
within the core of a compact star, it could provide 
a region within which rotational vortices are pinned 
and hence a locus for the origin of (some) pulsar 
glitches. Or, the presence of crystalline quark mat- 
ter within neutron stars could be ruled out if it pre- 
dicts glitch phenomenology in qualitative disagree- 
ment with that observed. 

There are two key microphysical properties of 
crystalline quark matter that must be estimated be- 
fore glitch phenomenology can be addressed. The 
first is the pinning force. Estimating this will require 
analyzing how the CubeX and 2Cube45z respond 
when rotated. We expect vortices to form, and ex- 
pect the vortices to be pinned at the intersections of 
the nodal planes at which condensates vanish. Ana- 
lyzing the vortices in three-flavor crystalline phases 
will be nontrivial. One complication is that because 
baryon number current can be carried by gradients 
in the phase of either the {us) crystalline conden- 
sate or the {ud) condensate or both, and the most 
favorable vortex or vortices that form upon rotating 
the CubeX and 2Cube45z phases will have to be de- 
termined. Another complication arises because the 
vortex core size, 1/ A, is only a factor of three to four 
smaller than the lattice spacing a. This means that 
the vortices cannot be thought of as pinned by an 
unchanged crystal; the vortices themselves will qual- 
itatively deform the crystalline condensate in their 
vicinity. 

The second microphysical quantity that is re- 
quired is the shear modulus of the crystal. After all, 
if vortices are well-pinned but the crystalline con- 
densate can easily deform under shear stress, the 
vortices will be able to move regardless of the pin- 
ning force. Glitches occur if vortices are pinned and 
immobile while the spinning pulsar's angular veloc- 
ity slows over years, with the glitch being triggered 



by the catastrophic unpinning and motion of long- 
immobile vortices. In order to immobilize vortices, 
and hence make glitches a possibility, both the pin- 
ning force and the shear modulus must be sufficient. 
The shear modulus can be related to the coefficients 
in the low energy effective th eory that describes the 
phonon modes of the crystal [33 . Isil l53 | . This effec- 
tive theory has been analyzed, with its coefficients 
calculated, for the two-flavor crystalline color super- 
conductor with face-centered cubic symmetry |5^. 
Extending this analysis to three-flavor crystalline 
color superconducting phases with the 2Cube45z 
and CubeX crystal structures is a priority for future 
work. 

Now that we have two well-motivated candi- 
dates for the favored crystal structure of the three- 
flavor crystalline color superconducting phase of cold 
quark matter, favorable over a very wide range of 
intermediate densities, the challenge becomes cal- 
culating the shear modulus and the pinning force 
exerted on rotational vortices in these structures. 
These arc the prerequisites to determining whether 
observations of pulsar glitches can be used to rule 
out (or in) the presence of quark matter in the crys- 
talline color superconducting phase within compact 
stars. 
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APPENDIX A: NEUTRALITY OF 
SOLUTIONS WITH Az = A3 

In Section VI. A, we gave a general analysis of 
the free energy ri(A2,A3). We showed that if we 
write (A2, A3) as \/2(Ar cos^, A^ sinf?) the free en- 
ergy takes the form H73() . and therefore has extrema 
only at 6* = 7r/4 (namely A2 = A3 = A^) or 
6 = 0, 7r/2 (namely a two flavor crystalline phase 
with only one A/ nonzero). As we have explained 
in Section II. D, in the strict Ginzburg-Landau limit 
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in which Aj/S^ any solution (A2, A3) is neu- 
tral. (The argument is that choosing = M'^/{Afj,) 
as in neutral unpaired quark matter suffices since, 
unlike BCS superconductivity, crystalline color su- 
perconductivity does not require any modification 
of the unpaired Fermi momenta prior to pairing and 
since in the Ginzburg-Landau limit the modifica- 
tions to number densities due to the pairing itself 
vanishes.) In this Appendix, we take a small step 
away from the strict Ginzburg-Landau limit. We as- 
sume that Aj. is small, but do not work in the limit 
in which it vanishes. We then show that the only 
solutions with = Mg/(Afj,) and, consequently, 
(5/i2 ~ = Sfj, ^ M^/(8/z) which are electrically 
neutral are those with A2 = A3 = A^. The two- 
fiavor crystalline phases with only one A/ nonzero 
are not neutral in three-flavor quark matter. 

The result of this Appendix allows us to neglect 
solutions which have only one A/ nonzero. This is 
fortunate, because there are many two-flavor crystal 
structures for which the sextic coefficient 7 is neg- 
ative, meaning that to sextic order the Ginzburg- 
Landau potential r2(A2, A3) often has runaway di- 
rections along the A2 and A3 axes [s^. Further- 
more, if the coefficient multiplying sin^ 9 in l\T6\\ is 
negative, for example if /3 and 7 are both negative 
while P32 and 7322 are both positive as is the case for 
both the CubeX and the 2Cube45z crystal structures 
on which we focus, then the extrcmum of 0(A2, A3) 
that we find with A2 = A3 appears to be a local 
maximum with respect to variation of 9 away from 
7r/4 while keeping fixed. We show in this Ap- 
pendix that upon fixing jUg = / (4/x) any solution 
with A2 ^ A3 is not neutral. For this reason, all 
these complications can be neglected, and we are 
correct to focus only on solutions with A2 = A3. 

The more formal way to proceed would be to 
define an ri„outiai(A2, A3), obtained by varying /ig 
(and ^3 and too) at a given value of the A's 
in order to obtain neutrality, and then finding A2 
and A3 that minimize r2noutrai(A2, A3). We have 
done a partial version of this investigation in a few 
cases and have found that, as expected, r^ncutrai docs 
have a minimum with /ig very close to /(4/i) and 
A2 very close to A3. A full exploration in this 
vein requires evaluating the Ginzburg-Landau co- 
efficients without assuming (5/i2 « Sji^ and, more 
challenging, requires reformulating our analysis to 
include nonzero and fig- We have not attempted 
the latter, and it is in this sense that our prelimi- 
nary investigation referred to above was "partial" . 



We leave this to future work, and turn now to the 
promised derivation of the neutrality of solutions 
with A2 = A3 and fj.e = M^/{4fj.). 

We shall only consider crystal structures for which 
{€[2} and {q^} are exchange symmetric, as this is 
the symmetry that allows the free energy to have 
extrema along the line A2 = A3. (Recall that by ex- 
change symmetric we mean that there is a sequence 
of rigid rotations and reflections which when applied 
to all the vectors in {q2} and {q^^} together has the 
effect of exchanging {(12} and {q2}-) Because we 
wish to evaluate dVl/dfie at = -^-'^^/(4/^)i we must 
restore /ig to our expression for the free energy Jl, 
rather than setting it to Af^/(4/i) from the begin- 
ning. Recall from 128|l that ficrystaiunc is the sum of 
the free energy for unpaired quark matter, which we 



know satisfies 9f2u 



cd/9/ig = 0atAtg = M2/(4Ai), 



and ri(A2, A3). Upon restoring the /Xg-dependence, 
the latter is given by 

r!(/ig,A2,A3) = 



Pa{6^l2)Al + Pa{5^l:i)Al 



-^32AlAl 



i i ^72 A2 + J^^^^Al + 7233 (^M2, Sti3)AlAl 



7322 ((5/^2, (^Ais) A3 A2 



(Al) 



where Sfi2 and Sfi^ can no longer be taken to be 
equal, as they are given by 



i5/^3 



2 

Ml 
4/i 



2 



which in particular means that 

9l5/Lt3 dSfl2 



dUe 



(A2) 



(A3) 



Because {£[2} and {0(3} are exchange symmetric, 
l32 = = P and 72 = 73 = 7- Because 
6^2 7^ <5/i3, however, the coefficients 7322 and 7233 
arc not equal and, furthermore, their {6(12, ^fJ-a)- 
dependence cannot be factored out as in H57|) or 
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(|7T|l . The cocfEcient 7322 depends on 81^2 and 
(5^3 through its dependence on JCudusus- 7322 = 

(^/^) X/ ^ud usu s (Qs; Q3! Q2 ' *^2i Q2 1 ^2 )• ^udusus IS 

given in (|67|l . Note that its dependence 

on i5/X2 and (5/i3 comes via q2 = ^rj6^2^2 

and = 5^*3 in addition to the ex- 

phcit dependence visible in (|67|l . Similarly 



7233 = (3/ 2)x;^«s«d«d(q2,q2>q3>qi>q3'q3) '^'^ere 

ICusudud has the same form as (|67|l except that (5/i2 
and 6fi3 are interchanged. Using the definitions H57(l 
and l(7T|l . one can confirm that (|A1|I reduces to (|7^ 
if we take 6^X2 = and hence 7322 = 7233- 

We now differentiate £7 given in I)A1|I with respect 
to /ie, noting ljA3p . obtaining 



917 _ dA2 

d^e 9A2 d/ie 
1 



9A3 rf/Xe 
1 



25^2^^!^ 2(5^2^M3 



1 ^ 97233(^/^2, ^Ms) 
3 



pda{5jiz)_^2^ _ p da{S^i2) ^2 



/332a2a: 
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97233 (<5Ai2,i5/i3) 



9(5/13 



9(5 /i2 



A^A^ 



d/i2 ^ ' 5/i| ^ 
J^A^ 

_^ 1 / dj322{6fi2,Sn3) 

3 \ dSfi3 9(5/i2 



97322(^^2,(5/^3 



(A4) 



We shall only evaluate d^l/dfj,e at values of A2 and A3 which are solutions to the gap equations 917/9 A2 = 
and 9S7/9A3 = 0, meaning that the first two terms in (|A4|) vanish. Furthermore, we shall only evaluate 
dfl/dfie at /te = Mg/(4/i), where 6^2 = ^/^3 = Sfi, and at solutions for which A2 = A3. Under these 
circumstances, the terms involving a, /32, /332 and 72 vanish and ljA4|l becomes 



on 

9/ie 



37r2 



97233 97: 



233 



97: 



322 



97 



322 



9(5/<3 9(5/X2 9(5/i3 dS^2 



A" ■ 

mill 



(A5) 



We argue that this vanishes as follows. Consider 
a particular term that contributes to 8^322/ d6fi2, 

dK.udususici3,ql^,qi,q,2,q.i,q.^)/d6n2- This is a 
complicated integral of a function which depends 
on the unit momentum vectors (cts, q2, qi, q2 , qi) 
and on dfi2 and dfis. From rotational invariance, we 
know that the value of the integral can depend on 
the relative orientation of the unit momentum vec- 
tors and on 5/Z2 and S^s but must be independent of 
common rotations of all the unit vectors. Now, all 
the crystal structures that we consider are exchange 
symmetric, meaning that for every quintuple of unit 
momentum vectors, (q-j, q2, qi, ^2 , q2) with the 
first chosen from {q3} and the last four chosen from 

{qa} there exists a quintuple (q2, qs, ql, ^3 ' qi) 
with the first chosen from {q2} and the last four 
chosen from {^3} such that the unit vectors in 
each of these two quintuples have the same relative 
orientation among themselves. This means that for 



every term 9K;„._(q|, qg, q^, qi, q^, qi)/95/.2. 

occurring in 8^322 /dSfi2, there is a correspond- 

ing term 9^_.„.(q^, q^, qf , qf, q{, qg)/9,5/Z3 

occurring in 97233/9(5 /is such that 

9/C„s«ciud(q2,q2,q3,q3>q3'q3)/9^M3 is related 

to 9/c„dus«s(q3,q3,q2,q2,q2,q2)/9^/^2 by the 

interchange of dfi2 and 6(13. Consequently, for 
Sfi2 = Sfi3 the two contributions cancel pair by 
pair when we evaluate 0^322/95^2 ~ 97233 /95/Z3 
or 97322 /9(5/^3 — 97233 /95/i2. In this way, the 
right hand side of l|A5|l vanishes, as we set out 
to show. We conclude that solutions to the gap 
equations with A2 = A3 and /ie = M|/(4/i) 
meaning (5/t2 = (5/^3 arc neutral. 

It is easy to see that the cancellations required 
in the proof of neutrality do not occur for solutions 
with A2 7^ A3. For example, following a derivation 
analogous to that above, we find that a solution with 
A2 ~ and only A3 nonzero is neutral with /ig = 
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M2/(4/x) only if 

a condition which has no reason to be satisfied. The 
study of sohitions with A2 7^ A3 therefore requires 
that they be constructed from the beginning with 
(5/i2 ^ (5/23 and with fixed by the neutrahty con- 
dition. We leave this to future work, focussing in 
Section VI on solutions with A2 = A3 which, we 
have proved here, are neutral. 

APPENDIX B: TRANSLATING {ua) 
RELATIVE TO {ud) DOES NOT AVOID 
REPULSION 

We have seen in Section VI that crystal structures 
in which a vector from {q2} and a vector from {q3} 
make a 180° angle are strongly disfavored, with infi- 
nite quartic and sextic Ginzburg-Landau coefficients 
(5^2 and 7322- Suppose we consider a structure like 
that in which {q2} and {q3} are coincident cubes, a 
disastrous choice. The way that we have improved 
upon this disastrous choice in Section VI. D is to ro- 
tate one cube relative to the other. Indeed, if we 
choose a 45° rotation about an axis perpendicular 
to a face of the cube, we obtain the 2Cube45z struc- 
ture which is one of the two crystal structures that 

I 



where we have worked only to cubic order. Using 
Q2 ~ Q2 + I2 = "12 conclude that the phase fac- 
tor in front of the J3i3i(q2, qii 121 12) term is sim- 
ply exp[2is(q|) • n]. In addition, we saw that for 
^3 ~ Q3 + ^2 = I2 to hold we need to have qg = q§ 
and q2 = qf- This makes the phase factor in front 



we find to be most favorable. In this Appendix, we 
ask whether wc can instead avoid the infinite free en- 
ergy cost of antipodal pairs by translating the {ud) 
condensate relative to the {us) condensate in posi- 
tion space, rather than rotating it. We find that 
the answer is no, and furthermore show that the 
Ginzburg-Landau free energy O that we have evalu- 
ated does not change if the {ud) condensate is trans- 
lated relative to the {us) condensate. 

Corresponding to each {q^} in momentum space 
we get a function A/(r) in position space which 
varies as A7(r) ~ X]q» e^*''^ To analyze the effects 
of translating A2 (r) relative to A3 (r) , it is helpful to 
restore the notation of ipiljl with A(qf ) representing 
the gap parameter corresponding to the momentum 
component q°. A2(r) or A3(r) can then be written 
as 

A,(r)=^A(q?)e2*'i°--. (Bl) 

Translating A2 (r) in the n direction by a distance s 
corresponds to the transformation A2(r) A2(r — 
sn) which multiplies each A(qf ) in the sum in I|B1(I 
by a different phase factor exp[— 2isq2 • n]. This 
is not just an (irrelevant) overall phase multiplying 
A2 (r) because it depends on the momentum compo- 
nent. The gap equation for the A2 components, as 
in (|40|l . is now given by 



(B2) 



of >/i2i3(q3, qs, q2^ ^2) s-lso exp[2is(q2) ■ n]. We con- 
clude that (up to cubic order) the gap equation for 
each A(q2) simply picks up an overall phase. The 
same is true for the gap equation for each A(q3). We 
therefore conclude that the free energy is unchanged 
up to quartic order when A2 (r) is translated relative 



A*(q-)e2-q^A^3l(q^,q^) 



^ A*(q^)A(q^)A*(q^)e2-('i^-q^+'i^)-J3i3i(q^,q^,q^,q5)5q._ 



\ A*(q^)A(q^)A*(q^)e2-'i^"Ji2l3(q^,q^,q^,q^)5, 



41 



to A3(r). This guarantees that such a translation 
cannot alleviate the large (3^2 arising from antipo- 
dal (or near antipodal) pairs of momenta occurring 
in {q2} and {Qs}- This argument can easily be ex- 



tended to include the sextic terms in the free energy; 
they too are unchanged when A2 (r) is translated rel- 
ative to A3(r). 
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